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CALCULATION OF SKIN-STIFFENER INTERFACE 
STRESSES IN STIFFENED COMPOSITE PANELS 


(ABSTRACT) 

A method for computing the skin-stiffener interface stresses in stiffened composite panels 
is developed. Both geometrically linear and nonlinear analyses are considered. Particular 
attention is given to the flange termination region where stresses are expected to exhibit un- 
bounded characteristics. The method is based on a finite-element analysis and an elasticity 
solution. The finite-element analysis is standard, while the elasticity solution is based on an 
eigenvalue expansion of the stress functions. The eigenvalue expansion is assumed to be 
valid in the local flange termination region and is coupled with the finite-element analysis 
using collocation of stresses on the local region boundaries. In the first part of the investi- 
gation the accuracy and convergence of the local elasticity solution are assessed using a 
geometrically linear analysis. It is found that the finite-element/local elasticity solution scheme 
produces a very accurate interface stress representation in the local flange termination re- 
gion. The use of 10 to 15 eigenvalues, in the eigenvalue expansion series, and 100 collocation 
points results in a converged local elasticity solution. In the second part of the investigation, 
the local elasticity solution is extended to include geometric nonlinearities. Using this analysis 
procedure, the influence of geometric nonlinearities on skin-stiffener interface stresses is 
evaluated. It is found that in flexible stiffened skin structures, which exhibit out-of-plane de- 
formation on the order of 2 to 4 times the skin thickness, inclusion of geometrically nonlinear 
effects in the calculation of interface stresses is very important. Thus, the use of a geomet- 
rically linear analysis, rather than a nonlinear analysis, can lead to considerable error in the 
computation of the interface stresses. Finally, using the analytical tool developed in this in- 
vestigation, it is possible to study the influence of stiffener parameters on the state of interface 


stresses. 
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1. Introduction 


In recent years there has been a dramatic increase in the use of composite materials in 
aircraft structures. Compared to metallic structures, the cost benefits, performance increase, 
and weight reductions which can be realized are substantial. As with metallic structures, one 
of the primary components in aircraft structures is the skin-stiffener combination. In traditional 
metallic structures the stiffener is attached to the skin by rivets. In severe loading conditions, 
such as those that occur during postbuckling, the rivets provide a site for yielding of the metal. 
As a result, failure in the form of local yielding may occur but the failure does not necessarily 
cripple the structure as a whole. On the other hand, composite materials are brittle and do 
not yield. Holes and other geometric discontinuities are sites for high stresses that ultimately 
cause the failure of the structure as a whole. Therefore, the use of rivets as the method of 
stiffener attachment is less attractive with composites than it has been with metals. A more 
common method of stiffener-to-skin attachment in composite structures is by an adhesive 
secondary bonding or by the cocuring of the stiffener and skin. 

Because of the differences in material properties and the lack of riveting, the failure char- 
acteristics observed with bonded or cocured stiffened composite skins are quite different from 
the failures encountered in riveted metallic structures. Failure in stiffened composite skins 
occurs in a much more catastrophic manner, being initiated by skin-stiffener separation [1J. 
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Further, it has been shown [2] that stiffened composite panels can fail prematurely, below the 
design load, due to skin-stiffener separation. One such mode of failure has been observed to 
initiate at the flange termination region. Figure 1 illustrates an example of this type of failure. 
Because of observations such as shown in Figure 1, it is believed that the stresses in the 
flange termination region are high. High skin-stiffener interface stresses can be attributed to 
a number of factors. First, there is the structural incompatibility associated with the deforma- 
tion of the skin and deformation of the stiffener when the skin-stiffener combination is sub- 
jected to applied loads. This incompatibility is particularly acute in the postbuckling state. In 
addition, if the stiffener is bonded to the skin, the flange termination region leads to a ge- 
ometric discontinuity in the structure. This region tends to serve as an area of increased 
stress. Finally, in stiffened composite skins the problem is compounded by an additional 
complication. In order to gain full advantage of the tailoribility of composite materials, gener- 
ally the stiffener and the skin are constructed of different material layups. Such material dis- 
continuity at the interface can lead to significant skin-stiffener interface stresses. It is the 
primary purpose of the present study to develop a method by which the stresses in this region 
can be accurately determined, and to use the method to investigate the influence of various 
stiffener parameters on the stresses in this region. Though the flange termination region is 
of primary concern, stresses at all locations along the interface are computed. 

To predict the skin-stiffener interface stresses, different levels of analytic complexity can 
be used. However, it is generally acknowledged that the analytical model should incorporate 
a few key features. First, the model should accurately represent the geometric and material 
discontinuities associated with the flange termination. As with free edge stresses, these ge- 
ometric and material discontinuities can cause the stresses in the flange termination region 
to be unbounded. When present, these unbounded stresses can be responsible for the initi- 
ation of skin-stiffener separation and should be accounted for. Secondly, since stiffened panels 
are most commonly designed to operate at the postbuckling range, the analysis must incor- 
porate geometrically nonlinear effects. Third, for the analytic tool to be useful in the design 
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Figure 1. Skin-Stiffener Separation Initiation at the Flange Termination Region. 


process, it must be sensitive enough to various design parameters, e.g., stiffener geometry 
and skin/stiffener material architecture. In addition to these three key factors, there is the is- 
sue of computational efficiency or cost. The computational cost should not be so large as to 
prohibit its use as a tool for parametric studies or for the design of a stiffened structure. And 
finally, the analysis should be such that it can be integrated into any proposed computational 
schemes or testbeds [3,4]. 

Considering the severity of the problem, there have been few studies of skin-stiffener 
interface stresses. One of the more notable ones is discussed in [5]. In this investigation the 
skin and flange were treated as separate orthotropic plate elements. These elements were 
held together via interface forces, the forces being taken as unknowns. The stiffness of the 
other stiffener elements, such as the web and the cap, were treated as extensional and rota- 
tional springs. The solution to the problem was formulated using the principle of virtual work 
and the theorem of minimum potential energy. Although in the postbuckling range the skin 
will actually experience moderate to large rotations, the analysis assumed only geometric 
linear deformation theory. Other studies [6-9] have considered similar problems, namely the 
adhesively bonded lap joints. In these studies the two adherents were treated as plates under 
cylindrical bending and/or inplane loads. In [6-8] the adhesive was modeled as shear spring 
only or tension-shear spring combination, and in [9] it was treated as an elastic layer in which 
the stresses did not vary across the adhesive thickness. In all of these studies a solution for 
the variation of shearing and peeling stresses in the adhesive was obtained. In [9] the derived 
solution was also applied to a simplified stiffened composite plate geometry. These investi- 
gations looked at simplified geometries and loading conditions. Hence, although these studies 
are important in furthering the understanding of adhesively bonded components, they cannot 
be applied in their present form to the study of stiffened composite aircraft structures. This is 
due to the fact that such structures tend to have complex geometries and loading conditions. 
In addition, the above studies did not address the issue of geometric nonlinearities. 
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Due to the problems in the previous analytical models discussed above, a more detailed 
analysis was undertaken. As stated above, and reiterated here, the following key require- 
ments were imposed: a) The model should accurately represent the state of stress near the 
point of geometric and material discontinuity; b) The model should be applicable to the ge- 
ometric nonlinear range; and c) The model should be sensitive in its stress prediction to var- 
ious design parameters, such as stiffener geometry and stiffener/skin material architecture. 
To meet these requirements, an analysis is developed in which the stress predictions are 
based on a generalized plane deformation elasticity solution in combination with standard 
finite-element calculations. The elasticity solution is valid in the localized region near the ter- 
mination of the stiffener flange. The finite-element calculations are valid for the structure as 
a whole, except near the flange termination region. The elasticity solution uses an eigenvalue 
expansion of the stress function to predict the stresses. The expansion is applicable in the 
flange termination region and satisfies exactly the boundary conditions there. The eigenvalue 
expansion is known to within arbitrary, but unknown, coefficients which are associated with 
each eigenvalue. The stresses from the finite-element solution and collocation scheme are 
used for determining the constants and thus uniquely determining the stresses in the localized 
region. 

In the next chapters various aspects of the investigation will be discussed. In chapter 2 the 
development of the analytical method and its verification are delineated. The derivations in 
this chapter are for geometrically linear analysis only. In chapter 3 it is shown how the meth- 
odology developed in chapter 2 may be extended to include geometric nonlinearities (i.e., 
moderate rotations). This is followed by a discussion of the application of this extension to the 
study of interface stresses in stiffened composite plates. In chapter 4, numerical results are 
presented. These results highlight the effect of including geometric nonlinearities in the 
analysis on skin-stiffener interface stresses. In addition, the effect of various stiffener ge- 
ometric and material parameters are evaluated. The study ends with some concluding re- 
marks and recommendations for future research. 
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2. Analytical Method Development 


2.1 Analysis Overview and Relevant Literature 


To facilitate the following discussion, a skin-stiffener cross-section being studied is shown 
in Figure 2. Note the location of the coordinate system in the figure and the nomenclature 
associated with the cross-section. Several coordinate systems will be used hereafter and it 
will be important to differentiate between them. Interest focuses on the computation of inter- 
face stresses at the skin-stiffener interface, along the line y=0. The type of stiffener shown, 
a blade stiffener, is only to serve as an example. Other stiffener types such as hat, I, and J 
can be studied with the type of analysis being developed. It should be noted that the flange 
can terminate at various angles, a, (see Figure 2), relative to the skin. The shaded area, 
shown in the exploded view, is referred to as the local region. This is the region where the 
elasticity solution is valid. The region outside of the shaded area, referred to as the global 
region, is the region where the finite-element solution is valid. It will be assumed that 
Figure 2 represents a cross-section of a stiffened panel that is long enough that the stress 
state does not vary in the stiffener direction in the region of interest. That is, in the nomen- 
clature of the figure, the stresses do not vary with the z coordinate. 
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Figure 2. Skin-Stiffener Cross-Section Geometry. 
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A linear analysis of the local region, as depicted in Figure 3, consists of finding the solution 
to the equations of elasticity valid in a region enclosed by an arbitrary corner opening created 
by two dissimilar orthotropic materials bonded along y = 0, and subjected to known boundary 
conditions. The literature regarding this problem is extensive and no attempt is made to cover 
all of it, instead some of the more important studies directly relevant to our investigation will 
be discussed. Williams [11] was the first to investigate the solution to an arbitrary corner in 
an isotropic plate in the form of an eigenfunction expansion. His solution was used extensively 
by Gross et al [12-16] in the numerical computations of stress intensity factors of various 
fracture toughness specimens. Carpenter [17] used Muskhelishvili's [18] complex potential 
approach to solve the same problem. He applied his solution to the investigation of interface 
stresses in lap joints. Williams, in a later paper [19], extended his solution to wedges 
composed of two dissimilar isotropic materials. Bogy [20] solved the same problem using 
Mellin transforms and obtained similar results. The extension of the above problem to wedges 
composed of two dissimilar anisotropic materials was first conducted by Wang and Choi [21]. 
In their work they used the complex eigenfunction expansion of the Lekhnitskii [22] stress 
functions. Later they applied their solution to laminate free-edge problems [23] and composite 
lap joints [24]. Other investigators solved the same problem using slightly different ap- 
proaches. Delale [25] transformed Lekhnitskii stress functions into polar coordinates and used 
a polar eigenfunction expansion, in real variables, to solve the problem. Ting and Chou [26] 
and Bogy [27] used the Green and Zerna complex function representation to solve the same 
problems. 

In the present study the solution method outlined by Lekhnitskii [22] and the specific 
eigenfunction expansion solution proposed by Wang and Choi [21] is employed. Inherent in 
this solution method is the need to determine a set of coefficients, one coefficient being as- 
sociated with each eigenfunction, from the boundary conditions. Among the many methods 
available by which the unknown coefficients may be determined, a collocation technique was 
chosen here. The method involves the use of n boundary conditions, of known magnitude, 
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Figure 3. Elasticity 


matched along prescribed boundaries with the same boundary quantities written in terms of 
the unknown coefficients. This leads to a system of n simultaneous equations from which the 
n unknown coefficients can be determined. The above procedure is sometime referred to as 
the point-matching method. Since there is no rigorous proof for convergence of the collocation 
method, it is typically studied by increasing the number of terms in the assumed function and 
the number of boundary collocation points. To increase accuracy and decrease dependency 
in the manner in which the positions of the collocation points are chosen, it is common prac- 
tice to use more boundary collocation points than unknown coefficients [28,29]. This approach 
is known as the overdetermined collocation procedure. Since it leads to an overdetermined 
set of equations by which the unknown coefficients are evaluated, the coefficients are deter- 
mined in a least-squares sense. Once it is clear that the number of terms in the assumed 
function and the number and location of the known boundary points have little influence on the 
numerical results, it is assumed that the procedure has converged. The converged functions 
are then assumed to be close to a true representation of the exact solution within this region. 

A number of investigators have made extensive use of the above described collocation 
procedure. Gross et al [12-16] used boundary collocation in the determination of K, for various 
edge crack specimen geometries. In their investigation the first coefficient {which is related 
to K,) of the Williams stress function, %, was determined by collocating boundary vales of x 
and (n being the normal to the boundary) along prescribed boundaries. Carpenter [28] 
applied boundary collocation in the determination of various fracture parameters. The un- 
known coefficients of the truncated function expansion were evaluated from the collocation 
of all three stress components ( o x , a y and t^) obtained from finite-element analysis. In an- 
other paper [30] Carpenter investigated accuracy issues related to the boundary collocation 
of stresses and/or displacements. Wang and Choi [23] used collocation to study the laminate 
free-edge problem. 

Other procedures which involve the use of eigenfunction expansion of the stress function 
in a localized region include the reciprocal work contour integral (RWCI) method, and the use 
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of special singular elements. The RWCI is based on Betti's theorem of elastic bodies. Two 
sets of forces (S, and S', ) and displacements (u, and u', ) of the same direction but not the 
same magnitude acting at selected points along the body boundary are in reciprocal equilib- 
rium. Use of the theorem provides a scheme by which the boundary value problem discussed 
previously may be solved. In addition, the RWCI leads to a path-independent procedure. The 
above technique have been used in the computation of K, and/or K„ by Carpenter [17] and 
Sinclair et al (31J. In the singular or hybrid element formulation, a special element is devel- 
oped which comprises a portion of the localized region where geometric and material dis- 
continuities occur. The displacements and stresses within the element boundary are 
governed by the exact elasticity solution. The surrounding standard elements of the mesh are 
then connected at n nodal points along the boundary of the special element. The coefficients 
of the truncated stress function are then determined so as to render continuity (or compat- 
ibility) of the nodal displacements at the special element's boundary in an exact or approxi- 
mate manner. There are numerous studies which have utilized the above concept. Some of 
the more relevant ones include work by Wang and Yuan [32] and Jones and Callinan [33], 

Clearly, the issue here is determining the conditions on the boundary of the local region. 
In general, it is not important how the boundary conditions are obtained, as long as the in- 
formation is accurate. For complex geometries, such as the skin-stiffener cross-section, the 
only reasonable method to obtain boundary conditions is with finite-element analysis. Thus the 
method here will utilized a finite-element analysis of the cross-section to provide stress in- 
formation on the boundary of the localized region. Furthermore, the finite-element analysis of 
the cross-section itself will be coupled to a finite-element analysis of the entire plate. The total 
analysis will be of the form of a structure-substructure-local analysis. The remainder of this 
chapter will be devoted to the local analysis. The governing equations of elasticity, the 
eigenvalue expansion solution, and the application of the collocation procedure will be pre- 
sented. In addition, the accuracy and convergence of the method are discussed by application 
to specific problems. 
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2.2 Material model 


The analysis here is implemented on the laminate rather than a lamina level. That is, the 
skin and stiffener are treated as having homogeneous integrated material properties. The 
reason behind this is as follows: It is felt that the interaction between the skin and the stiffener 
are controlled more by the overall stiffnesses of the skin and the stiffener than by the stiffness 
of the individual lamina at the skin-stiffener interface. For this reason the integrated material 
properties are used in the present analysis. Furthermore, only symmetric balance laminates 
are considered. In obtaining integrated material properties, the lamina principal material co- 
ordinates, denoted as the 1-2-3 coordinates, correspond to the transverse, thickness, and fiber 
directions, respectively. The x-y-z coordinates correspond to directions transverse, normal, 
and colinear to the stiffener, respectively (see Figure 2). The fiber angle, <p , measures the 
angle between the 3 and z axis, a positive rotation corresponding to rotation of the Tiber from 
the z axis toward the x axis. 

The well known laminate constitutive relations in the 1-2-3 system are written symbolically 
as, 


S 


[ 1 ] 


e, and cr, being the strain and stress vectors in the 1-2-3 system and S being the compliance 
matrix in the same system. The transformation of stress and strain from the 1-2-3 system to 
the x-y-z system leads to, 

o., = T, c x , [2. a] 

£i = T 2 e x , [2.b] 
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where T, and T* are the transformation matrices for stress and strain respectively. For a dis- 
cussion of these transformations see Appendix A. Substitution of eqs. 2 into eq. 1 leads to the 
lamina constitutive relations in the x-y-z system, namely; 


= S 0 y 


where, 


[S] = VCSIT, = 


S 11 

S12 

S -13 

0 

S 15 

0 

eg 
! CO 

s 22 

S23 

0 

S 25 

0 

S31 

S 32 

S33 

0 

S35 

0 

0 

0 

0 

S44 

0 

§46 

S 51 

S52 

S53 

0 

S55 

0 

0 

0 

0 

^64 

0 

S 66 


and Sj, = Sj,. In inverted form eq. 3. a is written as, 

o x — C e x 


[3. a] 


[3.b] 


[4] 


The integrated laminate properties are obtained by smearing the individual lamina properties 
throughout the thickness of the laminate. This is achieved by defining an average stress 
through the laminate thickness, h, i.e., 

= Y J-v h S x dy . [5] 

Substitution of eq. 4 into eq. 5 leads to the laminate constitutive relations: 
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[6. a] 


where the overbar is dropped from the stresses for convenience. The laminate stiffness 
components Ay are given by, 


A 


U 



c ij ( yk 


Vk-i ) 


[6.b] 


n being the number of laminae and y k and y k _-, being defined as the through-the-thickness lo- 
cations of the laminae interfaces. In inverted form; 

e x = [a] o x . [7] 

It should be noted that although each lamina is considered to be anisotropic in the x-y-z sys- 
tem, the laminate constitutive law is that of an orthotropic material, i.e., the smeared laminate 
properties are orthotropic. 


2.3 Elasticity Solution 


As mentioned at the outset, the analysis will be developed for the linear case and then 
extended to the geometrically nonlinear analysis case. What follows is the linear analysis 
development. 
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2.3.1 Development of the Stress Functions 


Consider an arbitrary semi-infinite corner composed of two dissimilar orthotropic materials 
bonded along y=0 (or 0 = 0), as depicted in Figure 3. If the body in question obeys the fol- 
lowing restrictions: a) The dimension in the z direction is much larger than the cross-sectional 
dimensions; and b) The external loads on the lateral surface do not vary with z, it is possible 
that the stresses, and hence the strains, are independent of the z-coordinate. Such a condi- 
tion is referred to as a generalized plane deformation. The ends of the body may be subjected 
to axial force, P„ twist M w , and moments, both about the x and y axes (i.e., M x and M y ). If such 
end loads are present, the state of plane deformation will exist at some distance from the ends 
in a manner consistent with St. Vernant's principle. Under the above conditions the stress 
equilibrium equations become, 


do x 

-f 

di xy 

dx 


dy 

8x xy 

+ 


dx 


dy 

3x xz 

4- 

^ T yz _ 

dx 


dy 


- 0 , 


= 0 , 


= 0 


[8.a] 


[8.b] 


[8.c] 


The strain-displacement relations are given by, 


e y " 


e z = 


du 

dx ’ 

3 

il 

dw , dv 
dy dz ’ 

[9.a,b] 

dv 

dy 

Yxz = 

dw , _du_ 
dx dz 

[9.c,d] 

dw 
dz ’ 

Yxy = 

du + _dv_ 
dy dx 

C9.e,f] 
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where the strains are a function of x and y only. The constitutive relations of eq. 7 for a 
homogenous orthotropic body under consideration can be written in full as, 



The general expressions for the displacement functions are obtained by a series of inte- 
grations and differentiations of eqs. 8, 9, and 10. The step-by-step details are given in [22]. 
In general, the constitutive relations are written in terms of the displacements using eqs. 9. 
The integration of three of these equations, (keeping in mind that a x is independent of z) and 
the satisfaction of the remaining three equations leads to the displacement functions in the 
general form, 


u = _ Bl - a ? 3 2 2 - b 4 y z + U(x,y) + co 2 z - co 3 y + u 0 , 

B 3 

V = - - 2 g 33 - z 2 + B 4 x z + V(x,y) + <o 3 x - z + v 0 , 


[11. a] 
[11. b] 


w = ( B, x + B 2 y + B 3 ) a M z + W(x,y) 

+ ©! y - © 2 X + w o • [lie] 

where Bi, (i = 1, 2, 3,4 ), are arbitrary constants of integration, co, , (i = 1,2,3) are rigid 
body rotations, and u 0 , v 0 and w 0 are rigid body translations. The unknown functions U(x,y), 
V(x,y) and W(x,y) must satisfy the following conditions: 
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dll _ 


= P^a x + p 12 a y + a 13 (B.|X + B 2 y + B 3 ) , [I2.a] 


— P 21 a x + P 22 a 23 (®1 x ®2 y ®3 ) 


dU + dV 
dy dx 


[I2.c] 


- P**« + B 4 y , 


[12. d] 


dW _ 


P 44 ” t yz ®4 X < 


[12. e] 


where P,, are the reduced stiffness coefficients and are given by, 


Pij “ a ij 


a l3 a j3 


, i.j - 1.2, 4, 5, 6 


The compatibility equations are satisfied identically for the above displacement field since it 
is derived from the strain-displacement relations. In addition it can be shown [22] that, 

o z - B : x + B 2 y + B 3 - -a^r(a 13 Cx + a 23°y) • ^14] 


The stresses which satisfy alt the aforementioned assumptions can be derived from two stress 
functions, F(x,y) and *F(x,y) . If the stresses are written in terms of these functions as, 


_ _ d 2 F _ d 2 F 

* x " 1? ' ° y dx 2 ’ Xy 


[I5.d,e] 
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then the stress equilibrium eqs. 8.a through 8.c are satisfied identically. The equations gov- 
erning F(x,y) and 4'(x,y) are obtained by substitution of the above stress relations into eqs. 
12.a through 12.e and the elimination of U, V, and W by differentiation. For an orthotropic 
material these equations become, 


P22 -77 + <2Pl2 + P66)-|tT + Pll T7 - 0 , 

<3x dx 3y <3y 4 


Cl6.a] 


- 

3x 2 dy 2 


= - 2B„ 


[I6.b] 


The decoupling of the two equations is a distinct characteristic of orthotropic materials (i.e., 
the equations are not uncoupled for anisotropic material). The equation governing F is ho- 
mogeneous while the equation for 4* involves a particular solution. The solution for F and the 
homogeneous equation for 4* have the form (22J, 

F = F ( x + py ) , [I7.a] 


4* = 4* ( x + uy ) , 


[I7.b] 


where p and u are parameters to be determined. For the local region near the vertex of the 
bimaterial corner the solutions for F(x,y) and 4'(x,y) are approximated in (21) as, 


F(Z) 


7 X+2 

c 

(X + 1) (X + 2) 


[I8.a] 


4*(Z) 


7 ( 5 + 1 ) 

D -7 

(5+ 1) 


[18. b] 


where, 


Z = x + py , [19. a] 

Z = x + uy . [19 b[ 
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and C and D are arbitrary constants. The substitution of F(x,y) and *F (x,y) into eqs. 16. a and 
16.b leads to, 


C X (X- 1) Z X ' 2 [p 22 + ( 2 P 12 + P 66 ) p 2 + PhM 4 3 - 0 , [20.a] 

D8Z 5_1 (P44 + P 55 « 2 ) = 0 . [20.b] 


The first equation is satisfied under the following conditions: 

a) p has one of four unique values given by the characteristic equation 

P22 + ( 2 P12 + p66 ) I* 2 Pi 1 fi 4 = 0 . [21] 


Such roots do exist and they are always complex or imaginary (for detailed discussion see ref. 
22). Considering these four values of p , F(x,y) is given by, 


F (e) (x,y) 


4 

1 

k = 1 


7 X+ 2 


(X + 2)(X, + 1) 


[ 22 ] 


b) p is arbitrary and, 


X = 0 , 1 


This leads to 

F (a) (x,y) = b 1 x 3 + b 2 x 2 y + b 3 x y 2 + b 4 y 3 

4- b 5 x 2 + b 6 x y + b 7 y 2 . [23] 

The superscripts (e) and (a) designates the eigenvalue expansion and auxiliary solutions, 
respectively. 


The second of eq. 20 is satisfied under the conditions: 
a) 0 has one of two unique values given by the characteristic equation, 
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[24.a] 


044 + u 055 - 0 . 


This leads to, 


u 1,2 


= ± i 


P44 

055 ’ 


[24.b] 


from which 'P(x.y) becomes, 


¥ (e) (x,y) 


2 z* + 1 

I D k — - 

k = i * (5+1) 


[25] 


b) u has an arbitrary value and, 


8 = 0 , 


leading to, 

^ a> (x,y) = b e x + b 9 y . [26] 

Finally the particular solution of eq 16. b is taken as, 

y (p) (x,y) = b 10 x 2 + b 1iy 2 . [27.a] 

In the above the C k and D k are arbitrary complex constants and b,, i = 1 9 are arbitrary real 

constants. The constants b 10 and b„ are not completely arbitrary, namely 

044 bio + P55 bn = - 2B 4 . [27. b] 

In addition, X and 8 are unknown parameters at this point. The total solution for the two stress 
functions is then written as the sum of the component solutions, i.e., 

F(x,y) = F (e) (x,y) + F (a) (x,y) , [28.a] 
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my) 


[28. b] 


= T (e) (x,y) + ¥ !a) (x,y) + ¥ p) (x,y) . 


2.3.2 Expressions for the Stresses and Displacements 


The Cartesian components of stress are derived from eqs. 15.a through 15.e as, 



<*x = 

S C kl x 2 k Z X k + G< a) , 
k = 1 

[29. a] 


Gy = 

i C k z£ 4- 0< a) , 

k = 1 y 

[29. b] 


T xy = 

£ ^k ^k + T xy » 

k = 1 y 

[29.c] 


T yz = 

— y d 7^ + <r(®) 

Z, L> k ^k + T yz 

k = 1 

[29. d] 


T xz “ 

2 D k o k Z 5 k + tg t 
k = 1 

[29. e] 


where the auxiliary stresses are given by: 



a (a) = 

2 b 3 x + 6 b 4 y + 2 b 7 , 

[30.a] 


a (a) = 

CJy 

6 b 1 x 4- 2 b 2 y 4- 2 b 5 , 

[30.b] 


11 

ss- 

k-* 

- 2b 2 x - 2 b 3 y - b 6 , 

[30.c] 


t w = 

T yz 

- b 8 - 2 b 10 x , 

[30.d] 

• 

•(•) = 
T xz 

b 9 + 2b„y . 

[30.e] 
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i 

It turns out that it is more convenient to impose conditions on the boundary of the localized 

i 

region using the polar cylindrical stresses. Using the cylindrical coordinates r - 0 - 

z (Fig- 


ure 3), the stress components become, 




cr r = I Ck^fmpk - n) 2 (m + 

k = 1 

np k ) x + a< a) , 

[31. a] 

4 

CT 0 = S c/(m + np k ) (X+2) 

k = 1 

+ 4 a) . 

[31. b] 

4 

4 X 

*r9 = E C k r (mp k - n)(m + 

k = 1 

n|./ +11 + T r0> • 

[31. c] 


t z9 = - I D k r s (m + nu k ) (5 + 1) + tg , 

k - 1 

[31. d] 

4 

2 c 

T rz = I D k r (mu k - n)(m + 
k = 1 

n u k ) 8 + t[ a) . 

[31.e] 

4 

where, 




o< a) = m 2 a< a) + 2mntJ + 

n 2 o^ a) , 

[32.a] 


4 a) = n 2 o<f> - 2mnx (a) + 

mX' 1 , 

[32. b] 

4 

T< a 0 > = n m (Oy 3 * - + (m 2 

- n 2 )t<J> . 

[32.c] 


4* e > = - niff + m,W , 


[32. d] 

4 

4z = . 


[32.e] 

4 

with m = cos 0 and n = sin 0 . 




The general displacement functions u, v and 
unknown functions U(x,y), V(x,y) and W(x,y) can 

w are given by eqs. 11. a through 11 
now be determined by substitution 

!.c. The 
of the 

4 
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stresses eqs. 29. a through 29. e into eqs. 12. a through 12.e. Integration of of eqs. 12. a, 12.b and 
12. e leads to, 


4 2^.+ 1 

U(x,y) = £ C k p k — ■ + U (a) (x,y) + U 0 , 

k = 1 (*■ + 1) 


V(x,y) = £ C k q k 


7 X+1 

^k 


k = l K (X + 1) 


+ V (a) (x,y) + V 0 , 


2 2^ + ^ 

W(x,y) = £ D k r k — + W (a) {x,y) + W 0 . 

k = i (5 + 1) 


where, 


Pk = PiiPk + Pi 2 • 


^k “ Pl2 Pk 


+ hi 


Pk 


r k = - 


J 44 


u k ' 


and, 


U'*> (x,y) = 3„ (b, x + 6 b„ y + 2 b 7 ) + p 12 ( 3 b, x + 2 b* y + 2 b 5 ) 


+ -^-(B 1 x + 2B 2 y + 2B 3 ) x + g(y) , 


V'*>(x,y) = p 1 ,(2b J x + 3b 4 y + 2b 7 ) + p„(6b 1 x + 2b J y + 2b ! ) 

+ x + B 2 y + 2B 3 ) y + f(x) , 

W (a) (x,y) = - [ P 44 (b 8 + 2b 10 x) + B 4 x]y + h(x) . 
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[33. a] 
[33. b] 
[33.c] 

[34. a] 
[34.b] 

[34.c] 

[35. a] 

[35. b] 
[35. c] 
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In addition, U and V must satisfy eqs. 12.c and W eq. 12.d. This leads to 


g(y) - C P 12 ) ^3 + 3 p 2 2 &i ~ ~~ y ~ ®i 3 y 2 + u 0 , [35. d] 

f(x) = [ (Pee — P 12 ) + 3 P^ b 4 — 3 x 2 + Pe6 ^6 x + ■ [35. e3 

h(x) = P55 bg x + W 0 . [35.f] 

The constants U 0 , V„ and W 0 may be dropped since they represents rigid body motion, terms 
which were already included in the general formulation of the displacement functions (see eq. 
11 ). 

2.3.3 Application of the Boundary and Interface Conditions 
2.3.3.a The Eigenvalue Solution 


An examination of Figures 2 and 3 illustrate the conditions that must be applied to the sol- 
utions to have the solutions satisfy the conditions of the skin-stiffener interface problems. 
Specifically referring to Figure 3, the surface represented by 0 = and 0 = — n are gener- 
ally free of any tractions. If pressure-loaded panels are being considered, these surfaces could 
be exposed to the normal pressure traction. However, the magnitude of this traction relative 
to the magnitude of the stresses generated within the material is negligible and can be con- 
sidered zero. Hence, one condition on the analysis is that the surfaces at 0 = a, and 
0 = - n are traction free. In addition along the line 0 = 0, the stiffener and skin are joined. It 
is the intent of the joining to provide a condition of no slippage along this line, i.e., the dis- 
placements are continuous across the interface. Finally, from stress equilibrium arguments, 
the stresses , and , are continuous across the interface. These conditions provide the 
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necessary equations for determining some of the constants, and hence the characteristic of, 
F(Z) and TfZ). 


The conditions of traction-free boundaries are represented by (see Figure 3), 

4 1) 0'.a 1 ) = 4e( r - a i) = 4 1 e( r - a i) = 0 , [36.a] 

4 2) ( r . ~ n) = 4e’ (r, - rt) = (r, - n) = 0 , [36.b] 

where the superscripts 1 and 2 designate material 1 (flange) and material 2 (skin) respectively. 
For a perfect skin-to-fiange bond the interface conditions along 0 = 0 require, 


aj^r.O) = 4 2) (r,0) , 

[37.a] 

-1 - — - 
CD-^ 

o 

II 

o 

[37.b] 

tS(r.O) - T ze> (r, 0) , 

[37,c] 

u (1) (x,0) = u (2) (x, 0) , 

[37.d] 

v (1) (x,0) = v <2) (x, 0) , 

[37.e] 

w (1) (x, 0) = w (2) (x,0) . 

[37. f] 


Finally, at the ends of the cross-section (see Figure 3) the following integral conditions are 
required to be satisfied, 


J J A T xz dx d v = 0 ■ 

J Ja V 

dx dy =0 , 

[38. a, b] 

JJ A a z dxdy = P z , 


dx dy = M x , 

[38.c,d] 

J j A a z x dx dy = M y , J J A 

(Ty Z X - 

*xz y) dx dy = M w . 

[38.e,f] 

! 
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It should be noted that as far as the above end conditions are concerned, it is possible to 
impose kinematic boundary conditions rather than force conditions, or it is possible to impose 
a mixture of the two. For example, rather than impose the area integral of a t over the body's 
ends to equal P x , it can be required that w = e 0 z, e 0 being an applied axial strain. 

The substitution of stresses, eqs. 31. a through 32. e, and displacements, eqs. 33. a through 
35.f, into the traction-free boundary conditions, eqs. 36. a and 36. b, and stress and displace- 
ment continuity conditions, eqs. 37.a through 37.f, places certain conditions on the stress 
functions, F(x,y) and ^(x.y). However, following the substitution of stresses and displacements 
into eqs. 36. a through 37.f and the application of variable separation to the resulting ex- 
pressions, it is evident that the conditions on the eigenvalue expansion part of the solution 
separate from the conditions on the auxiliary and particular parts of the solution. Therefore, 
the conditions on the eigenvalue expansion part of the solution (i.e., F (,l (x,y) and y ( *’(x,y)) are 
treated separately from the conditions related to the particular and auxiliary part of the total 
solution (i.e., F (,) (x,y), y w (x,y) and \|/W(x,y) ). The imposition of traction-free and traction and 
displacement continuity conditions on F (#) (x,y) and \|/ w (x,y) leads to the eigenvalue problem 
associated with 8 and X. Hence these parameters are determined uniquely for each problem. 
The conditions on F w (x,y), y (,) (x,y), and \|/ (p) (x,y) for the specific skin-stiffener configuration are 
discussed later, while the conditions on the eigenvalue expansion are discussed next. 

The traction-free boundary conditions associated with 'J'(x.y) are given by eqs. 36 as, 

4V(r.«i) = 4eV. ~ n) = 0 , [39.a] 

and the traction and displacement continuity conditions by eqs. 37 as, 

4e(r.O) = 4eV,0) , [39.b] 

w (1) (x,0) = w (2) (x,0) . [39. c] 

Equations 39. a through 39. c lead to the following relations associated with y^x.y): 
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[40.a] 


2 D[. 1) r S (cos a 1 + u{^ sin a 1 ) (5 + 1> = 0 , 

k = 1 


£ D< 2 >r 5 (-1)< 8 + 1 > = 0 , 


[40.b] 


k = 1 


(1) r 8 M >(8+1) _ 


2 Dk r (1) 
k = 1 


2 D< 2 >r 8 (1)< 8+1 > . 

k = 1 


[40.c] 


£ n (1) r * <5+1> 

«5, k k lsT^ 


l D«»r * (S+ ” 

k 5, Dk rk 7rnr ■ 


[40.d] 


where r k is defined by eq. 34. c. The above 4 simultaneous set of equations forms the 
eigenvalue problem for 8. This set of equations can be written symbolically as, 

[ Q (a, . 5) ] D = 0 . [41. a] 

For a nontrivial solution the values of 8 are given by, 

I Q (a, , 8) I = 0 . [41. b] 

The eigenvector, D , consists of the two eigenvectors, D (1 > and D m , related to materials 1 and 
2, respectively, i.e., 

— ► 

D = 

where both D (1) and D (a consist of two constants each. 


$<D 

;< 2 ) 


[41. c] 


If 8 is complex, solutions occur in complex conjugate pairs of the form, 

8 = y ± i P . [41. d] 

However, in order that w(x,y) be finite at the origin, 
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i 


t 

- 1 < y . 

[41. e] 

1 

The eigenvalue problem for 8 is described in greater detail in Appendix B. 



The boundary conditions associated with F (x,y) are given by eqs. 36 as 


4 

4 1 ’ O'. «i) = 4e ( f .«i) = 0 • 

[42.a] 


o^ 2) (r, - n) = (r, - n) = 0 . 

[42.b] 

4 

The traction and displacement continuity conditions associated with F(x,y) at the 

interface 


(0 = 0) are given by eqs. 37, 


4 

CT^fr.O) = 4 2) (r,0) , 

[43. a] 


4e (r.0) = 4e(r,0) , 

[43. b] 

4 

u (1, (x,0) = u (2) (x,0) , 

1 

| 

[43.c] 


1 

| v (1 * (x.0) = v< 2 > (x,0) . 

[43. d] 

4 

Equations 42. a through 43.d lead to the following 8 simultaneous equations associated with 


F'*»(x,y) : 

1 



£ C^ 1 * r* ( cos a 1 + p^. 1 ' sin a 1 ) (X+2) = 0 

k = 1 

[44. a] 

4 

2 cP^Oafc 1 * cos a, - sin a^){ cos a., + p^ sin a.,) tx+1 * = 0 

k = 1 

, [44.b] 

4 

0 

II 

oT 

+ 

1 

-C 

o 

^MII 

[44.c] 


£ c*V l tf><-i> |l+, > - o . 

k = 1 

[44.d] 

4 
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[44. e] 


£ 4 1) r x (1) (X+2) = £ 4 2) r x {1) (X+2> 

k=1 k=1 




k = l 


k = l 


[44.f] 


£ c<%(^ 


k = 1 


(X+1) = v c (2) (2) X (X+1) 

(X + 1) k 5i k Pk (X + 1) 


[44.g] 


y pH) n H) . ^ ^ — Y r( 2) n^ 2 ) * 

^ u k Qk TrjTTT " 4 Ik T 

= 1 (X + 1) k = 1 (X 


(X+1) 

(X + 1) 


[44.h] 


The above set of equations form the eigenvalue problem for X which is written symbolically 
as, 


[ A (a, , X) ] C = 0 . 


[45. a] 


For a nontrivial solution the values of X are given by, 


I A (a-, , X) I =0 


[45. b] 


The vector C is composed of C (1) and C (2 > associated with materials 1 and 2, respectively, 


C = 


id) 


i(2) 


[45. c] 


The eigenvectors C (1) and C ,2) consist of four constants each. If X is complex, solutions occur 
in complex conjugate pairs of the form 


X - r\ ± i£ . 


[45. d] 


However, in order for the displacements u(x,y) and v(x,y) to stay finite at the origin, 

- 1 < q . [45. e] 
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A detailed description of the eigenvalue problem for X is given in Appendix B. 


The eigenvalue problems associated with X and 5 do not occur in a standard form and hence 
require special procedures in order to determine the eigenvalues. Two methods were used in 
the present investigation, both methods involve the computation of the characteristic ex- 
pressions (in closed-form or numerically) which is equal to the determinant of the particular 
matrix of interest, i.e., Q of eq. 41. b or A of eq. 45. b. The first method, Muller (34,35], operates 
on a complex characteristic equation to find the roots. Once a root is found it is eliminated (or 
deflated) from the characteristic expression. This method is particularly well suited for com- 
plex root computation (see ref. 35). The second method is based on the secant technique for 
simultaneous nonlinear equations [36]. Here the real and imaginary parts of the determinant 
are treated as a set of two simultaneous equations, with the two unknowns being r\ and ^ or 
y and p . 

Using the relations between the stresses and the stress functions, eqs. 15.a through 15.e, 
the eigenfunction expansion of the stresses can be written [23], For the n* h real eigenvalue 8„ 


and X„, these stresses take the form, 

“ I c n [ £ Refc&rffzJ") + Imfc^+^VZ^) ] + o' a|> , [46.a] 

n=1 k=1 

oV = I C n [ £ Re ( Cpj*. z£" ) + Im (cU| k+2) z£ n ) ] + Oy a,) , [46.b] 

*xy = ~ I c n [ £ Re(cMz£") + ImtcfJk+j^Z^) ] + . [46.c] 

7 n = 1 k = 1 7 

= - I d n [Re(d*!Ui n ) + lm(dM")] + . [46.d] 

n = 1 7 

tS - - Z d n [Re(d<^u ( 1 l) Z^) + + t£'> . [46.e] 

n = 1 
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Whereas, for the n ,h complex eigenvalues 8„ and X„ the stresses take the following form: 



= 


xz 


0) JO 78 


= I { d n Re[ S dfjfeoji 

n = 1 k = 1 


z-] + 


d' n Im [ £ 

k = 1 


dW: 


] } + 


r ( a i) 

l xz 


[47.a] 


[47.b] 


[47.c] 


[47.d] 


[47.e] 


where, i = 1,2 corresponding to material 1 (flange) and material 2 (skin). In addition, 
c n ,c' n ,d„ and d'„ are unknown real coefficients, whereas, c*J£ and d|;£,i = 1,2, are known 
quantities of the normalized n^ eigenvector associated with materials 1 and 2 respectively. 
The reader is referred to Appendix B for a detailed discussion of the appropriate eigenfunction 
representation for a real and complex eigenvalue. Next the application of the boundary and 
interface conditions on the auxiliary and particular portions of the solution is discussed. 


2.3.3.b The Remaining Part of the Solution 


In previous sections we saw that the imposition of traction-free boundary and interface 
conditions on the eigenvalue expansion part of the solution for F(x,y) and 'P(x.y) led to the 
eigenvalue problems for X and 8, respectively. In this section these conditions are enforced 
on the other parts of these solutions (i.e. F w (x,y) , 4 #(,) and (x,y)). The imposition of the 
boundary and interface conditions are implemented in the context of the skin-stiffener geom- 
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etry as discussed in section 2.3.2. The traction and displacement continuity conditions are 
given by eqs. 37. These conditions require that, 

u m (x,0) = u ( 2 ) (x,0) , [48. a] 

v ( 1 ) (x,0) = v ( 2 ) (x,0) , [48. b] 

4y(x.O) = xg?(x,0) , [48. c] 

Tyz(X.O) = 'tyz(X.O) , [48. d] 

Oy^x.O) = <jy 2 ^(x.O) [48. e] 

Here we chose to use the Cartesian components of stress. Considering eqs. 11 . a through 11 . c 
and 29. a through 30.e the above relations become: 

- + uM (x , 0| + 4 <> z + „<<> = [49, aj 

h (2| J2I 

- S±J*L Z * + U W (X , 0 ) + c4» z + u l=l , 

R (1) a (1) 

- 2 2 33 z 2 + B 4 x z + v (ai) (x,0) + ©y'x - ©V’z + = [49. b] 

r(2) _(2) 

- ~ 2 — 2 2 + B^xz + V (aj) (x,0) + 4 2, x - <o^ 2) z + v< 2) , 

(BV’x +B< 1 >)a' 1 3 ) z + W (ai, (x,0) - 4 1l x + w' 1 * = [49.c] 

(B^> x + B^ 2) ) ag z + W (a ^(x,0) - 4 2) x + w< 2 > , 

- 2 4 1) x + b^ = - 2b ( 2 2) x + b^ 2) , [49.d] 
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[49.e] 


2 b\g x + bj, 1 * = 2b!, 2 0 ) x + bg 2) 

6 b^x + 2 bf* = 6b!j 2) x + 2b^ 2) . [49.f] 


Matching coefficients of the same power of x,y, and z we arrive at the following relations: 


E>!”a£> = BP'ag 

, i = 1,2,3 , 

[50,a] 

Bi” - 

CD 

[50.b] 

of" = coP' . 

i = 1,2 , 

[50.c] 

II 

*> 

i = 1,2,5,6,8,10 

[50.d] 

U (a,) (x,0) + u* 1 * = 

U (aj) (x,0) + u< 2) , 

[50.e] 

V (ai) (x,0) + co^x + v^ = 

V (a?) (x,0) + 4 2) x + v‘ 2 > . 

[50.f] 

W (a,) (x,0) + w' 1 * = 

W (aj) (x,0) + w« 2) , 

[50.g] 

where the superscripts 1 and 2 correspond to material 1 (flange) and material 2 (skin). Next 
the traction-free conditions are applied to the skin and flange free surfaces. These conditions 
for the skin, eq. 36. b, require that, 

<4 2) (r, - 7i) 

= 0 , 

[51. a] 

tr9 O'. “ *) 

= 0 , 

[51. b] 

4o ( r ' “ *) 

= 0 . 

[51.c] 


Consideration of eqs. 29 through 32 in conjunction with the above conditions leads to the fol- 
lowing relations: 
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i 


- 6b!| 2, r + 2b^ 2) =0 . 

[52. a] 

i 

2 bj 2) r - bg 2) = 0 . 

[52. b] 


b^ 2) - 2b< 2 >r = 0 . 

[52. c] 

l 

Equating coefficients of the same power of r yields 



b[ 2) = 0 , i = 1,2,5,6,8,10 . 

[53] 


Finally, considering eq. 50.d together with eq. 53 leads to 



bj 1 ’ = bj 2) = 0 , i = 1,2,5,6,8,10 . 

[54] 

1 

Based on the above results, the auxiliary stress components, eqs. 30, reduce to 



o^ a) = 2 b 3 x + 6 b 4 y + 2 b 7 , 

[55. a] 

« 

O 

f! 

b 

[55. b] 


t< a y > = - 2 b 3 y , 

[55. c] 

i 

o 

II 

H>s 

[55. d] 

% 

= b 9 + 2 b„ y . 

[55.e] 


Further information can be obtained by considering the traction-free conditions for the 


flange (material 1), These conditions are 


i 

oiVai) = 0 , 

[56. a] 


T r 0 * (r. a-|) = 0 , 

[56. b] 

t 
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i 


4e (r. a i) “ 0 


[56. c] 


The use of these conditions, together with eqs. 32 and 55, leads to the following relations: 


• 

6rn 2 [mt4 1) + nb£ l) ] + 2bp ) n 2 = 0 , 

[57.a] 


- 2rn[bp ) (2m 2 - n 2 ) + Sb^mn] - 2mnb| 1) =0 , 

[57.b] 

• 

- n bg ' - 2 b!jV r n 2 = 0 , 

[57.c] 


where, m = cos(a ,) , n = sinfa,) . Equating coefficients of the same power of r leads to 

• 

bp 1 = 0 , i = 3,4,7,9,11 . 

[58] 


Finally, considering eqs. 50.e through 50.g in conjunction with eqs. 27. b, 35, 54 and 58 leads 

• 

to 



0«> b «> = f bV> - f B® . 

[59.a] 

• 

o R (2)u(2) = a (1) R (1) - a (2) B (2) 

‘ Pi 1 d 7 a 13 b 3 a 13 °3 

[59. b] 


3 P fi'bf - *jpBf> - f b!?> , 

[59. c] 

• 

4 2) = o , 

[59. d] 

# 

9 d(2) 

b (2> - - 2B * 
$ 

[59.e] 


Applying relation 50.a produces 


• 

bf - KB® , b® - • bf - KB® . 

[60.a,b,c] 
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where, 


K 


J1)J 2) _ J2)J1) 

a 13 a 33 a 13 a 33 




[60.d] 


Finally, considering all the relations derived in this section, the auxiliary stresses for the skin 
and stiffener are as follows: 


Flange (material It 

oW = o< ai) = T< a ;> = T< a '> = t£> = 0 . [61. a] 

Skin (material 2) 

a*f 2> = 2 K (b!| 2) x + Bj, 2) y + B^ 2) ) , [61. b] 

= 0 , = 0 , [61.c,d] 

o r(2) 

>&’ - - 2KBfy . tg>> = jj-y . [ei.e.f] 

Pm 1 


The stresses in material 1 (flange) and material 2 (skin) are given by eqs. 46. a through 47.e. 
These stresses are given in terms of the unknown coefficients c„, c'„, d n , d'„, and 
B,. i = 1,2, 3, 4. For the specific skin-stiffener problem these coefficients are determined by a 
collocation procedure. This is discussed next. 
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2.4 Collocation Procedure 


The eigenfunction expansion elasticity solution is valid for a semi-infinite domain. The 
particular boundary value problem associated with the skin-stiffener geometry is solved by 
assuming that the solution is valid in the finite domain represented by the localized flange 
termination region. Furthermore, the solution is assumed to be represented by a finite number 
of terms in the series. The unknown coefficients c n and c'„ in eqs. 46.a through 46.c and 47.a 
through 47. c are determined using the boundary collocation technique. Referring to 
Figure 2, the boundary of the local region is bounded by contour ABCDEFA. In the actual 
structure analysed, boundaries AF, AB, BC, and ED are traction free, whereas boundaries CD 
and EF are subjected to both normal, o n , and tangential, x, tractions. Obviously in the nomen- 
clature of the problem, o n = <r x and x, = x w on CD and EF, and <j„ = a y and x, = x^ on BC and 
DE. By the development of the elasticity solution, the traction-free conditions on FA and AB 
are already satisfied. The collocation procedure is used to satisfy the traction-free conditions 
on boundaries BC and DE, and to match the normal and tangential tractions on boundaries 
CD and EF as determined by the global finite-element analysis. It should be noted that in the 
boundary collocation procedure, the normal and shear stresses, a„ and x, , along the contour 
BCDEF are written in terms of the unknown coefficients, c„ and c' n in a truncated eigenfunction 
expansion. These stresses are matched with the same stress components calculated by the 
finite-element analysis on the contour. Although it is possible to collocate other responses, 
such as three components of stress, the strains, the displacements, etc., it is felt that matching 
the normal and shear stresses on the boundary is the best choice. The main reason for this 
is the fact that enforcement of the force equilibrium conditions on the local region as a finite 
body involves only the normal and shear traction on the boundary of the body. 

For the specific problem here, the collocation at m points around boundary BCDEF leads 
to 2m simultaneous equations from which the 2n unknown coefficients are determined. The 
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use of more collocations points (i.e., m > n) leads to an overdetermined set of equations. 
Solution of these equations produces the 2n unknown coefficients c„ and c' n which satisfy all 
boundary conditions in a least-squares sense. If the original set of equations is represented 
by, 


S = A C , [62. a] 

then the least squares solution for C is [37], 

C = (A t A)“ 1 A t S . [62.b] 

In the above, A is a 2m x 2n known matrix (for which m > n ), S is a vector of length 2m 
consisting of a 2m known boundary stress quantities, and C is a vector of length 2n consisting 
of 2n unknown coefficients, c„ and c'„ , of the truncated eigenfunction. The elements of matrix 
A involve material properties and the coordinates of the collocation points. Once these coef- 
ficients are determined, the stresses in the localized region can be written. For a converged 
eigenfunction, these stresses are assumed to represent the true stress field in this region. In 
a later section convergence is studied by varying the number of terms (or eigenvalues) in the 
truncated function expansion and by varying the number of collocation points. 


2.5 Global Finite-Element Analysis 


To facilitate the application of the local elasticity solution to the skin-stiffener problem, a 
finite-element program was implemented which incorporated the generalized plane deforma- 
tion assumption. This finite-element formulation is consistent with the elasticity solution de- 
veloped earlier. The finite-element program was developed for two reasons. The primary 
reason was to provide boundary conditions for the eigenvalue expansion local elasticity sol- 
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ution. The second reason was to check the results of the eigenvalue expansion. As indicated 
in the introduction, the methodology was developed for geometrically linear problems, and 
then extended to geometrically nonlinear problems. Thus the finite-element program imple- 
mentation was very important in the first step of the solution strategy development. 

A body which conforms to the generalized plane deformation assumption may be analysed 
using a two-dimensional finite-element model. In this section a short description of the finite- 
element program (PE2D) is given. The PE2D program developed is based on the FEM2D 
finite-element program, [38], with the appropriate modifications to meet the need of the pres- 
ent investigation. For a more detailed discussion the reader is referred to Appendix C. 

The generalized plane deformation elasticity finite-element model is based on the dis- 
placement field given by eqs. 11. a through ll.c for homogeneous anisotropic bodies for which 
stresses do not vary along the generator (i.e., the z axis). This displacement field can be 
written in vectorial form as, 


u = u 0 + U 


[63] 


where; 


and, 


U T = {U(x,y) ; V(x,y) ; W(x,y)> , 

Uq = {u 0 (x,y,z) ; v 0 (x,y,z) ; w 0 (x,y,z)> , 


u 0 (x,y,z) = - - - 1 |^ z 2 - B 4 yz 

V 0 (X,y,Z) = _J2|33. Z 2 + b 4><z 
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w 0 (x,y,z) = (B^ + B 2 y + B 3 )a 33 z 


Here the terms associated with the rigid body rotation and translation are omitted. The un- 
known functions U(x,y), V(x,y) and W(x,y), are approximated according to the finite-element 
method. Further, it may be shown that the unknown constants Bj , (i = 1 ,2,3,4) are related 
to the body's kinematic end conditions, that is: 


B 3 - 


e o 

^33 * 


B 2 




[64] 


In the stiffened skin structural context, e 0 and k, are the axial extension and curvature in the 
z direction, and k„ is the twist curvature about the z axis. The coefficient B, is related to the 
inplane twist about the y-axis and is of no consequence in the structure considered here and 
is therefore set to zero. The strain vector, e*, is given by, 


= E + e r 


[65] 


where; 


5 


e t 


rdU 
1 dx 


-T 

o 


= {0 ; 


{ 8 X • I 

E z > Yyz 

> Txz 

* Txy } 

a 

>1 

. cW , 

<3W 

dU , 

dy 

’ dy ' 

dx 

dy 

0 ; (e 0 + k z 

y) ; - 

K xz 

2 X 1 

Kxz 
2 y 


Since only orthotropic media are being considered, the constitutive relation is as given before 
(see eq, 6. a) as 


o x = [a] e* x . 
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Considering eqs. 6.a and 65, it may be concluded that the problem of determining V(x,y) and 
U(x,y) decouples from the problem of determining W(x,y). In addition, since U,V and W are 
functions which depend on x and y alone, only a two-dimensional finite-element model is re- 
quired. The finite-element model is derived via the variation formulation in a standard manner. 
The conditions for this model consist of specified overall kinematic conditions (i.e., e„, k„ and 
k „ specified ) and force and/or kinematic conditions on the boundary in the x-y plane. It 
should be pointed out that because of the generalized plane deformation assumption, the 
specified overall kinematic conditions e 0 , k z and k« do not. vary along the the z-axis of the 
body. 


2.6 Verification of the Analytical Model 


In this section attention is given to a results relevant to the verification of the linear analysis 
model. In addition, the fidelity of the local-global elasticity-finite-element analysis is demon- 
strated. By fidelity is meant the accurate representation of stresses by the elasticity solution 
within the localized region and the smooth transition to the global region. Whereas the con- 
vergence of the finite-element method has been discussed in [39], there is no rigorous proof 
for the convergence of the collocation method and its application to the present problem. For 
the problem here the accuracy of such a procedure will depend on a number of factors such 
as: a) the accuracy of the boundary conditions; b) the number of eigenvalues used in the 
eigenfunction expansion; and c) the number of collocation data points. In the following section 
these issues will be addressed. 

In order to study convergence and accuracy of the current model, the elasticity results are 
compared with results from a finite-element analysis in which the mesh was refined twice in 
the localized region. The results are compared for a particular problem. The use of a fmite- 
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element analysis for comparison is due to the lack of another analysis of this particular 
problem. Two flange geometries are considered in this phase of the study, a 90° and a 45° 
flange termination angle, i.e., a, = 90° and a, = 45° in Figures 2 and 3. Two angles were con- 
sidered at this stage so as to make conclusions regarding model verification more general. 
The flange and the skin are constructed of an 8-ply quasi-isotropic ( ± 45/0/90), laminate, and 
an 8-ply orthotropic ( ± 45/90 2 ), laminate, respectively. Material properties are given in Ap- 
pendix D. The finite-element discretization for 90° and 45° flange angle skin-stiffener geom- 
etries are shown in Figures 4a and 5a, respectively, for what is referred to as the coarse mesh. 
Subsequent mesh refinements of the localized region are shown in Figures 4b, c and 5b, c. Due 
to the geometric and material symmetry, only one-half of the structure cross-section is mod- 
eled. The particular problem considered for the verification study is shown Figures 4.a and 5.a, 
namely the stiffened plate subjected to a pure bending moment M. The particular loading was 
chosen to illustrate the computational method because this loading produce peeling and 
shearing stresses at the skin-stiffener interface that are approximately the same order of 
magnitude. It should be mentioned that for this problem and are identically zero. 
Therefore, attention is focused on o x , a y , and t^. The exact nature of the skin-stiffener inter- 
face stresses near the flange terminus depends on the the value of ti,, the real part of the first 
eigenvalue (see eq. 45.d). These stresses are unbounded if -1 < 

In the collocation procedure, the normal, cr n , and tangential, t„ stresses are collocated 
along the closed contour ABCDEFA shown in Figures 4 and 5 by the heavy line. The stresses 
as computed by the finite-element analysis of the entire cross section are used to provide 
collocation data on the contour. In reality, collocation takes place only along boundaries BC, 
CD, DE, and EF, since the conditions of stress-free boundaries along AB and FA are satisfied 
exactly by the elasticity solution. Boundaries BC and DE were taken as stress-free faces of the 
skin and flange, respectively, whereas the stresses along the internal boundaries CD and EF 
are those determined by the finite-element analysis. These stresses on the boundary are 
calculated exactly within the finite-element context by postprocessing the finite-element dis- 
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(a) Coarse Mesh 



(b) Refined Mesh 



(c) Fine Mesh 


Figure 4. Loading and Finite-Element Discretizations of Skin-Stiffener 90 s Flange Termination 
Angle. 
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placement data. The collocation points are uniformly spaced along the contour. The bounda- 
ries CD and EF are placed at distance L from point A. This distance, shown in Figures 4 and 
5, was chosen according to the condition that a y and as computed by the finite-element 
analysis must be continuous across the skin-flange interface. That is, the stresses o y and 
computed by the elements on either side of the interface in the skin and flange must yield the 
same value to within 5%. The theory of elasticity stipulates that these two components of 
stress are exactly continuous across the interface. This 5% condition assures that the skin- 
stiffener interlaminar stresses computed by the finite-element satisfy this condition to with a 
small tolerance and thus are accurate in the global region. Based on this 5% condition, it was 
determined that for the coarse mesh, L should be about 1.5 times the combined thickness of 
the skin and flange. This value of L was kept constant for all subsequent mesh refinements 
in order to eliminate variation in the results which may depend on L. In the subsequent fig- 
ures the stresses are plotted as a function of the normalized distance, x = x/t„ from the flange 
termination vertex (point A in Figures 4 and 5), where t, denotes the skin thickness. In all fig- 
ures, whenever a x stress data is presented, values computed both in the skin and in the flange 
are given. This dual computation is done because the stress component a x is discontinuous 
at the interface and will have different values in the skin and flange. However, if interface 
stress data for o y and t^, as computed by the elasticity solution, are given, only one value is 
provided, since as stated above, these stress components are continuous across the interface. 
On the other hand, for the finite-element computation of these two stress components, both 
the stresses in the skin and the stresses in the flange are given. Plotting both the flange and 
skin finite-element components will illustrate the degree to which the continuity condition on 

these stresses is violated as the stress gradient becomes severe near the flange termination 

M t 

region. Finally, all stress components are normalized by — p where M is the applied mo- 
ment, -y is the distance from the neutral axis of the skin to the outer surface of the skin, and 
I is the moment of inertia of a section of skin of unit depth into the paper. This strength-of- 
materials view of normalizing the stresses is meant only as aid to not having to become in- 
volved in discussing the actual magnitude of the numbers. The normalized stresses are 
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denoted as a „, a y , and Throughout the study, the stresses, and infact, distances, have 
been normalized in a manner relevant to the particular problem. 

Figures 6 through 11 show the variation of the skin-stiffener interlaminar stresses, a x , a y , 
and as a function of distance along the interface for the 90° and 45° flange termination an- 
gles. Both stresses determined by the finite-element analysis and those calculated using the 
truncated eigenfunctions, eqs. 46 and 47, are shown. The elasticity results shown in the figures 
were generated using 15 eigenvalues and 100 collocation points. The eigenvalues for both 
geometries are given in Appendix D. Each figure illustrates calculated values of a particular 
stress component for each of the three different meshes (i.e., coarse, refined, and fine). The 
data point closest to the vertex for which finite-element stress data is plotted is at x = 0.025. 
It should be noted that x goes from 0 to 3 because the skin and flange are of equal thickness, 
t, , and the length of the local region, L, is 1.5 times the combined thickness of the flange and 
skin, 2t,. The scale of the vertical axes was taken to account for the value of the normalized 
stress at x = 0.025. 

A number of interesting observations can be made from these figures. In general, for both 
flange termination angles, at sufficient distances from the flange termination vertex, the 
finite-element and the elasticity solutions show excellent agreement for all three components 
of stress. The point at which this agreement can be categorized as being excellent moves 
closer and closer to the vertex with increases in the mesh refinement. As will be discussed 
shortly, there is very little difference between the elasticity solution resulting from collocating 
data from a coarse mesh and the elasticity solution resulting from collocating data from the 
fine mesh. In essence, for either flange termination angle, the elasticity solution does not 
change from one mesh to the next and hence the use of the coarse mesh is sufficient. What 
the figures are showing is that the elasticity solution coincides with the solution the finite- 
element analysis is appearing to converge to. It should also be noted that o y and computed 
by the finite-elements adjacent to the interface are discontinuous across the interface in some 
region near the vertex. The situation is worst for than for a y . This discontinuity relates to 
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FINITE ELEMENT, FLANGE 
FINITE ELEMENT, SKIN 
ELASTICITY, FLANGE 
ELASTICITY, SKIN 



Figure 6. Skin-Stiffener Stress g x for 90° Flange Termination Angle, Finite-Element and Elasticity 
Results Using IS Eigenvalues and 100 Collocation points. 




FINITE ELEMENT, FLANGE 
FINITE ELEMENT, SKIN 
ELASTICITY 



Figure 7. Skin-Stiffener Stress a y for 90° Flange Termination Angle, Finite-Element and Elasticity 
Results Using 15 Eigenvalues and 100 Collocation points. 









FINITE ELEMENT, FLANGE 
FINITE ELEMENT, SKIN 
ELASTICITY 



Figure 8. Skin-Stiffener Stress for 90° Flange Termination Angle, Finite-Element and Elasticity 
Results Using 15 Eigenvalues and 100 Collocation points. 








FINITE ELEMENT, FLANGE 
FINITE ELEMENT, SKIN 
ELASTICITY, FLANGE 



Figure 9. Skin-Stiffener Stress <7 X for 45* Flange Termination Angle, Finite-Element and Elasticity 
Results Using 15 Eigenvalues and 100 Collocation points. 
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FINITE ELEMENT. FLANGE 
FINITE ELEMENT. SKIN 
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Figure 11. Skin-Stiffener Stress z xy for 45° Flange Termination Angle, Finite-Element and Elasticity 
Results Using 15 Eigenvalues and 100 Collocation points. 






the inherent inaccuracy of the finite-element data near points of stress singularity. Although 
the region of inaccuracy shrinks with increases in mesh density near this point, the discrep- 
ancy between the two values of stress near the vertex is quite significant. This characteristic 
occurs for both the 90° and the 45° flange termination angles. 

To illustrate that the stresses as computed by the coarse mesh are sufficiently accurate to 
serve as boundary conditions for the elasticity solution, Figures 12 and 13 illustrate the vari- 
ation with x of the three stresses a x , o r and as generated by the elasticity solution for the 
two flange geometries. Each figure illustrates the three components of stress and each stress 
component is represented by three relations, each relation corresponding to the collocation 
of stress data from the different finite-element meshes. It is immediately evident from the fig- 
ures that for each stress component for both flange angles, the results of collocating stress 
data from the three different mesh densities coincide. This indicates that in order to obtain a 
good approximation of the state of stress in the localized region using the elasticity solution, 
only the coarse mesh need to be used. This is somewhat expected since at some distance 
away from the point of singularity the coarse mesh will yield a converged set of stresses. That 
is, the values of the stresses which were used for collocation on the boundaries CD and EF 
were approximately the same for all meshes. This supports the computation effectiveness of 
the current procedure, since only a coarse finite-element mesh is required to produce highly 
accurate stresses in the localized region. From the figures it is hard to distinguish between 
the magnitudes of the stresses at x = 0.025. Hence, the values of ct„, ct y , and at these 
points, for the three meshes, are given in Table 1. It is observed that the percent difference 
in these values is negligible, again emphasizing the power of this methodology. 

At this point it is appropriate to depart from the discussion of the accuracy of the method 
and illustrate some of the physical results that can be addressed from the analysis. Table 1 
indicates that in terms of the peel stress, the 45° flange termination angle is better than the 
90° flange termination angle. At x = 0.025, the normalized peeling stress is 0.290, while for the 
45° angle, the peeling stress is 0.192, 33% less. On the other hand, the shear stress is about 
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COARSE MESH 
REFINED MESH 



Figure 12. Skin-Stiffener Stresses for 90° Flange Termination Angle Computed From the Elasticity 
Solution using Three Mesh Refinements, 15 Eigenvalues, and 100 Collocation points. 







COARSE MESH 
REFINED MESH 
FINE MESH 



Figure 13. Skin-Stiffener Stresses for 45° Flange Termination Angle Computed From the Elasticity 
Solution using Three Mesh Refinements, 15 Eigenvalues, and 100 Collocation points. 
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the same for the two different flange geometries. On balance, the 45° angle would be better 
of the two. This coincides with intuition but the methodology developed here allows for the 
behavior of the two different angles to be quantified. 

Next, the issue of the convergence of the truncated eigenfunction expansion, as related to 
the number of eigenvalues and collocation data points, is addressed. Figures 14 and IS illus- 
trate the variation with x of o x , a y , and as computed by the truncated stress functions (eqs. 
46 and 47) for 90® and 45° flange termination angles. The multiple data on each plot corre- 
spond to stress calculations produced by 5, 10 and 15 eigenvalues in the truncated 
eigenfunctions. In these computations 100 collocation points from the coarse mesh were used. 
As can be seen in the figures, convergence of the truncated stress functions for o y and 
occurs between 10 and 15 eigenvalues. That is, the stress computation for 15 eigenvalues is 
bracketed between the computation for 5 eigenvalues and the computation for 10 eigenvalues. 
On the other hand, this convergence is not completely evident for a„. There seems to be a 
slight increase in this stress value with an increasing number of eigenvalues. There does not 
seem to be a bracketing effect with increases in the number of eigenvalues, as was observed 
for the other stress components. However, on a percentage basis, the increase of stress is 
minimal. Hence it is felt that the use of between 10 to 15 eigenvalues leads to convergence 
of the eigenfunctions in eqs. 46 and 47. Figures 16 and 17 illustrate data for the same stress 
components in which the number of collocation points were varied. In particular, the number 
of collocation points was doubled. For these computations 15 eigenvalues and a coarse mesh 
were used in the stress computations. From Figures 16 and 17 it is apparent that the number 
of collocation points has a minimal (or no) effect. This is true provided that the number of 
collocation points are approximately 3 to 4 times larger than the number of coefficients in the 
truncated eigenfunctions. It should be noted that the number of collocation points for the 90° 
flange termination angle, Figure 16, is different than the number of collocation points for the 
45° flange geometry, Figure 17, because of slight differences in geometry. 
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Figure 14. Elasticity Solution of Skin-Stiffener Stresses for 90° Flange Termination Angle as 
Function of the Number Eigenvalues Using 100 Collocation points. 
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Figure 15. Elasticity Solution of Skin-Stiffener Stresses for 45 s Flange Termination Angle as 
Function of the Number Eigenvalues Using 100 Collocation points. 
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Next a point of computational importance is addressed. In the course of the investigation 
it has been found that the stresses on contour ABCDEFA, more specifically the stresses on 
CD and EF, as computed by the finite-element analysis, are practically independent of the 
flange termination angle. Thus the stress data from the finite-element analysis of the cross- 
section with a flange termination angle of 90° can be used as collocation data for the elasticity 
analysis of a 45° flange termination angle, or any other angle for that matter. Hence, the 
finite-element analysis need to be done only once. This reduces the computation task, since 
the one finite-element analysis can be used to produce multiple local elasticity solutions for 
different flange termination angles. Figure 18 shows results of two elasticity solutions of the 
90° flange termination angle. In one solution the finite-element data from the analysis of a 
cross-section with a 90° termination angle and the eigenvalue expansion solution for the 90° 
flange angle are used. In the other solution the finite-element data from the analysis of the 
cross-section with a 45° termination angle and the eigenvalue expansion for the 90° angle are 
used. Of course the former solution is the one that is correct but it is clear from the figure that 
the latter solution is essentially identical to the former. Figure 19 shows the solution for 45° 
flange termination angle solved using both 90° and 45° finite-element analyses as the bases 
for the collocation procedure. Again it is obvious there is very little difference in the results, 
indicating that the eigenvalues and eigenvectors play a very important rule in determining the 
characteristic of the stress distribution in the local region. 

It is appropriate at this time to address an issue that is often raised in conjunction with the 
study of interface stresses. The point is raised now because it will put into context the results 
to be presented in the remainder of this document. The analysis developed to this point has 
assumed that the skin and the flange are joined along a line that represents a perfect bond 
with zero thickness. Depending on the fabrication process, the bond line may have nonzero 
thickness and the material within this thickness may have properties significantly different 
than either the skin or the flange materials. At issue is the fact that the interface between the 
flange and the skin may be an adhesive with nonzero thickness. Of concern is the influence 
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Figure 18. Skin-Stiffener Stresses for 90° Flange Termination Angle Using Boundary Data From 
the Finite-Element Analysis of a 90 s and 45* Flange Termination Angles. 










of the adhesive layer on the peeling and shearing stresses. In what follows, the interface 
stresses are computed using finite-element analyses for two adhesive layer thicknesses. 
These results are compared to results for the case for which a zero bond line thickness is 
assumed. The data presented provide important information as to whether the elasticity sol- 
ution is conservative or nonconservative in the calculation of skin-stiffener interface stresses 
when an adhesive layer is present. One may postulate that an elasticity solution in which the 
two adherents are assumed to have a perfect bond of zero thickness will lead to noncon- 
servative stress calculations. This is a consequence of the fact that the assumption of perfect 
bond reduces the flexibility of the interface, as compared to if the adhesive layer was present. 
Hence it is important to know how large an error is introduced as a consequence of the as- 
sumption made. To determine this effect, the same skin-stiffener geometry and loading con- 
ditions which are shown in Figure 4 are analysed. Two adhesive layer thicknesses are 

evaluated, - = 0.05 and y- = 0.10. The thicker adhesive layer represents the upper limit of 

•» ** 

the bond line thickness. The thinner layer thickness is the nominal bond line thickness. The 
results from these analyses are compared to stress data for y- = 0 layer thickness as com- 
puted by the finite-element program and as computed by the local elasticity analysis. The 
adhesive material properties are given in Appendix D. Mesh refinement in the localized region 
corresponds to the fine mesh of Figure 4.c. A coarser mesh is used outside the local region. 
However, element size in the coarse region is such that the adhesive elements' aspect ratio 
does not exceed 1:4. Figure 20 shows the character of the interface peeling stress o y in the 
localized region. Figure 21 displays similar information for the interface shearing stress t^. 
Each figure consists of two plots, one illustrating the stresses at the interface as computed 
by elements in the flange, and the other illustrating the same stresses but evaluated by ele- 
ments in the skin. The finite-element stress data is denoted by symbols, whereas, the local 
elasticity results are indicated as solid curves. As in previous figures the stresses and dis- 
tance along the interface have been normalized. 
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Figure 20. Skin-Stiffener Interface Peeling Stress, a , Throughout the Flange Termination Region 
as a Function of Adhesive Thickness. 










In general, the results of Figure 20 suggest that there is a slight increase in the magnitude 
of the peeling stress, a y , between x = 0.025 to x = 1 as the bond line thickness increases. 
However, the increase relative to the zero bond line results is quite small, something on the 
order of less than 5%. In general, the results of Figure 20 indicate that even in the present 
of adhesive layer, the peel stresses tend to become unbounded as the flange terminus is ap- 
proached. The unbounded nature of the peel stress with an adhesive is in qualitative agree- 
ment with the finding of [9], The conclusions regarding the behavior of the shear stress are 
not so clear. As has been shown in Figures 8 and 11, the elasticity solution and the finite- 
element results for zero thickness bond line are in disagreement near the flange termination 
point. Further, the finite-element shear stress results in the skin and in the flange at the 
interface near the flange terminus are markedly different from each other. This essentially 
indicates that the finite-element solution is not converged at this location. However, with the 
adhesive present, the finite-element stress calculation for the shear stress in the flange are 
not that different than the calculation for shear stress in the skin. Has the finite-element sol- 
ution converged and the shear stress is indeed finite and may be less than the zero bond line 
case ?. Reference [9] indicates that the shear stress is also quite large near the flange ter- 
mination point. This finding, coupled with what appears to be the poor shear response of the 
8-node finite-element for this class of problems would lead to the conclusion that the shear 
stresses are unbounded near the flange termination point, even in the present of an adhesive 
layer. Hence, it is felt that the peeling and shearing stresses computed by the elasticity sol- 
ution for the case of zero bond line thickness are representative of the behavior, even in the 
present of nonzero bond line thickness. 

Attention is now focused on extending the methodology to include geometric nonlinearities, 
and to study the influence of stiffener geometry and material properties on interface stresses. 
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3. Application of the Methodology to Stiffened 
Composite Plates 


3.1 Preliminary Remarks 


If the cross-section being studied is actually part of a plate which is subjected to a particular 
loading and boundary conditions, the cross-section experiences a loading more complicated 
than the simple bending moment used in studying accuracy and convergence. In addition, the 
loading is most certainly of such a magnitude that the plate experiences geometrically non- 
linear effects. Thus, to extend the present work to the analysis of stiffened plates with realistic 
loading and boundary conditions, the cross-section geometry being consider is assumed to 
be part of a pressure-loaded plate which is clamped on all four sides. The particular plate 
geometry and boundary condition were chosen since they represent the configuration studied 
experimentally in [10] and, as shown in Figure 22, the situation simulates a unit cell of stiffened 
aircraft skin structure. Although in actual aircraft structures the skin is generally subjected to 
both transverse pressure and inplane loading, only a pressure load is investigated. In that 
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Figure 22. Single Stiffener Plate Geometry. 


context it should be mentioned that the methodology can be extended to include more general 
loadings and stiffener geometries. In the present investigation a coordinate system different 
than the one commonly use by plate theory was chosen (see Figure 22). This was done in 
order to keep the coordinate system consistence with the elasticity solution developed earlier. 
Although, this may lead to some confusion, it was felt that the plate analysis (herein called the 
structural analysis) had only sideline importance in the present development, the elasticity 
analysis being the principal analysis. 

The analysis procedure to be discussed consists of three steps. First, a two-dimensional 
structural level analysis of the entire stiffened composite plate is conducted. Next, a typical 
cross-section in the central region of the stiffened plate structure is isolated and a three- 
dimensional finite-element analysis of the cross-section is performed. The boundary condi- 
tions for this three-dimensional model are provided by the structural level analysis. Hereafter, 
the finite-element analysis of the stiffener-skin cross-section is alluded to as the substructural 
analysis. In the previous chapter this type of analysis was referred to as global analysis. Here, 
it is felt that since the stiffened-skin cross-sectional analysis is a subset of the structural level 
analysis, this terminology is more appropriate. Finally, the stresses from the substructural 
analysis are coupled with boundary collocation scheme to produce a rigorous elasticity sol- 
ution in the flange termination region, i.e., the localized region. The scheme thus may be 
termed as a structure-substructure-local analysis. Both the structural and substructural ana- 
lyses were conducted using the commercially available Engineering Analysis Language (EAL) 
finite-element code [40]. Before proceeding to discuss the structure-substructure-local analy- 
sis model, the extension of the elasticity solution to account for geometric nonlinearities is 
considered. 
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3.2 Geometrically Nonlinear Elasticity Analysis 


Since stiffened composite aircraft skin structures are generally designed to operate in the 
postbuckling or geometrically nonlinear range, this characteristic must be considered. De- 
spite the geometrically nonlinear deformations, it can be assumed that the strains are small 
and hence linear elastic material response will prevail. In addition, for a thin flexible structure, 
such as a stiffened composite plate, it is customary to assume that the inplane deformations, 
u and w, are much smaller than the plate dimensions in x and z directions, whereas, the out- 
of-plane deformation, v, can be of the same order of magnitude as the plate thickness: The 
above assumptions lead to neglecting all nonlinear terms in the strain-displacement relations 
which contain the derivatives of u and w (" indicates inplane displacement of the plate 
midsurface). In addition, if it is assumed that a line perpendicular to the midsurface remains 
perpendicular and unstrained after deformation, then one arrives at the well-known Von 
Karman plate kinematic relations. 

The extension of the geometrically linear elasticity solution previously developed to the 
geometrically nonlinear range relies on (for the problem studied here) the assumption of small 
strains and small to moderate angular rotations, ©„ and co„ about the x and z axes. In addition, 
it is required that these rotations are independent of spatial location within the localized re- 
gion. Under these assumptions the stress equilibrium equations written in the deformed body 
coordinates are satisfied by the stress function relations, eqs. 15. a through 15.c, written in the 
same coordinate systems. To facilitate our understanding as to why such a statement can be 
made, a brief discussion of the general nonlinear stress equilibrium conditions follows. For a 
detailed discussion of the subject the reader is referred to Novozhilov, [41]. 

In the geometrically nonlinear range the equilibrium equations are derived by summing 
forces on an infinitesimal element in the deformed state. For the purposes of the analytical 
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development, it is sufficient to concentrate on the geometric representation of a body under- 
going small strains and small to moderate rotations. In the case of small straining, the 
infinitesimal element shape in the deformed state may be approximated by its original 
undeformed shape. In addition, the x-y-z coordinate system will be transformed into a rec- 
tangular cartesian system, x-y-z, in the deformed body, as shown by Figure 23. 

In the skin-stiffener nonlinear analysis context, the following assumptions are made: 

(a) The condition of small strains exists in the localized region and linear elastic material 
response is obeyed. 

(b) The condition of small to moderate rotations exists and these rotations are inde- 
pendent of spatial location within the localized region. 

Under the above considerations, the localized skin-stiffener region in its original shape 
undergoes rigid-body rotations as shown in exaggerated form in Figure 23. a. Since x-y-z', is a 
rectangular cartesian system, the equilibrium equations can be written in the deformed state 
by adding tildes over x,y and z. Hence, the elasticity solution developed previously is assumed 
to be valid, provided that it is applied in the deformed body state and that the angular rotations 
are spatially uniform within the region of its application. The procedure outlined previously 
for the determination of the coefficients, c n and c' n in the eigenfunction expansion, eqs. 46 and 
47, is then valid provided it is performed in the deformed body state. 


3.3 Structural-Substructural-Local Analysis Procedure. 


The structure-substructure-local analysis procedure is depicted by Figure 24. First, a 
structural analysis of the entire blade-stiffened plate (Figure 24.a) is conducted using plate 
elements. Next, a substructural analysis of an isolated region, where skin-stiffener interface 
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(a) structure model (b) substructure model (c) local model 



Figure 24. Structure • Substructure - Local Analysis Procedure. 


stresses are expected to be high, is performed (Figure 24.b). Finally, a local elasticity analysis 
of the flange termination region is performed (Figure 24.c). 

The structural finite-element model uses plane stress membrane and bending plate ele- 
ments. The substructure finite-element model employs three-dimensional brick elements. All 
of these elements are available in EAL. Using these elements, both linear and nonlinear an- 
alyses can be performed. The nonlinear solution algorithm used by EAL employes the full or 
the modified Newton-Raphson method. In the present investigation the full Newton method 
was used. Most of the elements in the EAL library (and all of the ones used here) are based 
on Pian's [42] hybrid element formulation derived according to the minimum complementary 
energy principle. In general terms, this formulation assumes a stress field (generally in 
polynomial form) which satisfies stress equilibrium in the element interior. The displacements 
along the element boundaries are expressed in terms of compatible generalized nodal dis- 
placements, in a manner identical to the one used in displacement formulation fmitfe-element 
development. Next, using the principle of minimum complimentary energy, the element 
stiffness matrix can be constructed. For a detailed discussion of this method see (42) and 
Appendix C. In the next sections each step of the analysis procedure is discussed. 


3.3.1 Structural Analysis 


3.3.1 .a Details of Finite-Element Analysis 


The stiffened plate is discretized using 4-node plate elements, referred to as the E43 ele- 
ment in EAL, as shown in Figure 25. Each element of this type has 5 degrees of freedom 
namely; 
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Figure 25. EAL Finite-Element Discretization of a Stiffened Plate Configuration. 




To verify of the structural analysis procedure and the accuracy of the EAL finite-element 
model, a stiffened plate configuration which was tested and analysed in [10] was reanalyzed 
using EAL, E43 elements. In [10] the STAGS finite-element model was used in the analysis of 
the plate. The stiffened-plate finite-element discretization used was similar to the one shown 
in Figure 25. The comparison between the EAL analysis, the STAGS analysis, and the exper- 
imental data is shown in Figure 26. In the figure the out-of-plane deflection, at the center-of- 
plate and in the skin away from the center of the plate, are plotted as a function of the applied 
pressure. Both the results from STAGS and from EAL are based on geometrically nonlinear 
analyses. The figure clearly show that the two finite-element methods give nearly identical 
results. For an explanation of the difference between the experimental data and the analytical 
results, the reader is referred to [10]. Results like the one displayed in Figure 26 provide 
confidence in the accuracy of the EAL structural analysis model used in the present investi- 
gation. 


3.3.1. b Details of Structural Analysis 


To be consistent with the constitutive law used in the elasticity solution, in the structural 
level analysis all extension and bending-twist coupling terms are set to zero. These terms 
include B 1g , B 2g , D ie , and D ze . For symmetric balanced laminates, such as the ones considered 
here, B ie = B^ = 0 for both the skin and the flange. On the other hand, neither D ie nor D 2e 
are zero. However, they tend to be small relative to D^, D 12 , D a , and D M . Given these material 
restrictions, and the geometry, and loading conditions considered, it was only necessary to 
analyse one-quarter of the plate using the two planes of symmetry (i.e., [x.y.O] and [O.y.z]). 

Before proceeding, an important point must be raised. The elasticity analysis of the flange 
termination region relies on the assumption that the stresses do not vary with the z coordi- 
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nate, i.e., the stiffener direction. The stresses in a stiffened plate which is clamped on all four 
sides clearly vary with the z coordinate, particular near the ends of the stiffener. However, 
toward the center of the plate there will be a region within which the stresses are quite uni- 
form with respect to the z coordinate. It is in this region that the elasticity analysis is valid. 
Consequently, it is in this region that the structural analysis will be refined to a substructural 
analysis, and the substructural analysis used to provide boundary data for the collocation 
scheme. It is fortunate that in many instances skin-stiffener separation occurs near the center 
part of the plate because the analysis methodology discussed here is indeed valid in this re- 
gion. 

Since attention will be focused on the central region of the plate, it is of interest to study 
the convergence of the structural level finite-element analysis in this area. Figure 27 illus- 
trates a coarse mesh and a mesh which was refined twice in the central region. These two 
meshes and one which was refined once, not shown in the figure, were used to study plate 
response, and in particular, convergence. In the mesh refinement, triangular elements were 
used as transition elements. The triangular element, E33, and the regular element, E43, are 
discussed in Appendix C. 

Figure 28 shows the convergence characteristics of the three components of displacement 
along x at z = b/3. The distance of z = b/3 is the distance from z = 0 to the second nodal line 
along x and b is the flange half-width. The second nodal line was selected over the first one 
since along the first nodal line, at z=0, both w and p x are zero and therefore cannot be used 
for convergence study. The displacements studied are the displacements at the skin 
midsurface. The figure shows the displacements for the three meshes used to study conver- 
gence, i.e., coarse , refined, and fine meshes. It is clear that the results for the refined and 
the fine meshes coincide, indicating convergence of the response along the line z = b/3. 
Figure 29 shows similar results for components of rotation, P x and P z along the line z = b/3 
Although the figures present displacements and rotations data at z = b/3, similar convergence 
characteristics were observed for the entire mesh refinement region. 
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(a) Coarse Mesh 



(b) Fine Mesh 

Figure 27. 1/4 Symmetry Stiffened Plate Diacretizatlon and Mesh Refinements. 
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with 10 psi Pressure (Data Plotted at z— b/3). 




COARSE MESH 
REFINED MESH 
FINE MESH 



Figure 29. Rotation Convergence Study in a Stiffened Composite Plate Nonlinear Analysis with 
10 psi Pressure (Data Plotted at z=b/3). 





One of the unique characteristics of stiffened composite plates is revealed in Figure 28. This 
is the phenomenon of pillowing. Pillowing is the term used to describe the fact that the skin 
at x/b = 3 deflects out-of-plane more than the center-of-plate, x/b = 0. This pillowing literally 
pulls the skin away from the flange. Pillowing is the reason the displacement of the skin in 
Figure 26 is greater than the displacement of the center-of-plate. 


3.3.2 Substructural Analysis 


The substructural analysis is conducted using 8-node hexahedrons and 6-node 
pentahedrons, solid elements which are available in EAL. Each node of these elements has 3 
degrees of freedom, namely u, v, and w. Both elements accept fully anisotropic material 
constitutive behavior. To be consistent with the elasticity solution developed, the material 
compliance matrix is taken to be orthotropic. A brief description of the two solid elements 
used is given in Appendix C. 

The substructural analysis is conducted on isolated area in the structure. As mentioned 
above, attention is given to the plate center location. The region for which the substructural 
analysis is carried out is indicated by the area outlined by the heavy line in Figure 30. In ad- 
dition, Figure 30 also shows schematic discretization of the substructure model. The model 
includes both sides of the symmetry plane and has a width of 4b/3. One side of the model 
coincides with nodal lines in the structure model at z = 2b/3. The model is discretized with 8 
elements in the z direction and the mesh refinement of the flange-skin cross-section corre- 
sponds to the coarse mesh used in section 2.6. As was done in the structural analysis, the 
web is modelled with plate elements except that now the web is attached to the top of the 
flange. The following displacements boundary conditions are imposed on the substructure 
model: 


3. Application of the Methodology to Stiffened Composite Plates 


85 



(a) Structure Model. 



(b) Substructure Model. 


Figure 30. Structural to Substructural Analysis Procedure and Schematic of Substructure 
Discretization. 
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Boundary 

Specified Displacements 

(x,y,-2b/3) 

u, v, and w 

(x,y, + 2b/3) 

u, v, and w 

(3b,y,z) 

u, v, and w 

(0,y,z) 

u = 0 

(x.y.O) 

w = 0 


where u, v, and w are determined from the structural analysis displacement and rotation data 
by applying Kirchhoff assumptions. That is: 


• 

u(x,y,z) 

= u 0 (x,z) - ~ 2 - y , 

[68. a] 


v(x,y,z) 

= v 0 (x,z) , 

[68. b] 

e 

w(x,y,z) 

3v n 

- w 0 (x.z) + - 5r y . 

[68. c] 


Here, u 0 , v 0 , and w 0 are the displacements of the skin midsurface as determined from the 
structural level analysis. On the -z boundary the displacement boundary conditions are eval- 
uated using the oddness and evenness of the displacement and rotation functions. These 
conditions are 

u 0 (x,z) = u 0 (x, - z) ; v 0 (x,z) = v 0 (x, - z) ; w 0 (x,z) = - w 0 (x, - z) , 


dv 0 (x,z) 

dx 


_ 


3v 0 (x, - z) 


dx 


dv 0 (x,z) 

dz 


dv 0 (x, ~ z) 

dz 


[69] 
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The boundary conditions determined by using eq. 68. b lead to unnatural through the thickness 
constraints at the boundaries of the substructural region. This restriction may be alleviated 
by specifying linear variation of v through the thickness, as was done in [44]. In the present 
study the local elasticity analysis was performed along the substructure symmetry plane 
(x,y,0]. Consequently, it was felt that the effect of through-the-thickness constraints at the 
boundaries, {z = ± 2b/3), would have little if any effect on the state of stresses of the location 
z=0. In addition to the displacement boundary conditions of eqs. 68, the substructure model 
was subjected to transverse pressure. Here the pressure load was applied in the form of a 
vacuum to the skin side. This was done in order to simulate the loading conditions used in 
110 ]. 

In a linear analysis, a set of displacements corresponding to a particular pressure level, 
as compute by the structural analysis, is imposed on the substructure model. The system of 
equations representing the state of the body is then solved once. Solution for any other 
pressure level may be obtained by simple ratio of the previous results. In a nonlinear analysis, 
a more complicated procedure is required. Specifically, the response at each pressure level 
is the result of a series of finite-element calculation. First, a converged set of displacements 
and rotations for a particular pressure level, say p 1 , is obtained using the structural analysis. 
This set of displacements and rotations are then used to compute the boundary conditions 
discussed above for the substructure model. Next, using these conditions, the set of equations 
representing the substructure is solved iteratively until convergence is reached for that pres- 
sure level. The converged solution represents the solution for one pressure level, p v The 
process must be repeated if another, say a higher, pressure level is required. Using the con- 
verged solution, the stresses can be computed at any point throughout the interior of the 
substructure using the stress polynomials which are part of the EAL finite-element formulation 
(see Appendix C). These stresses are then used in the collocation scheme. 

Before proceeding to discuss the local analysis procedure, it is instructive to compare the 
response predicted by the substructural analysis with the response of the structural analysis. 
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Probably the single most significant kinematic variable in the analysis of flexible stiffened 
plate structures is their out-of-plane deformation. Consequently, this variable was chosen as 
a measure of comparison between the two models. The results of such comparison are shown 
in Figure 31. In this figure the out-of-plane displacement v is plotted as a function x for both 
the structure and substructure models. The results correspond to nonlinear analysis using 10 
psi applied pressure. One portion of the figure corresponds to v(x,0) and the other portion to 
v(x,b/3), z = b/3 being half the distance between the symmetry plane and the boundary in the 
substructure model. The substructure extends to 3b in the x direction, as can be seen from the 
point of termination of the dashed line. The figure demonstrates the excellent agreement be- 
tween the two models. However, as would be expected, the substructure model appears to 
be slightly stiffer. One may improve on the agreement between the two models by adding 
additional elements in the z direction to the substructure model. This, off course, would in- 
crease the cost per run. The error observed with the present results are on the order of 1 to 
2 % and do not justify the increase in cost. 


3.3.3 Application of the Local Elasticity Analysis 


In practice there is only one difference between the geometrically linear and nonlinear local 
elasticity analyses. The difference between the two is that in the nonlinear analysis the 
collocation procedure must be performed in the deformed body configuration. In practice, this 
presents no problems. It should be pointed out that should collocation of displacements or a 
mixture of displacements and stresses been used, extension to the nonlinear case would not 
have been as simple. Given the discussion of section 3.2, the only requirement on the local 
geometrically nonlinear elasticity analysis which is different from the linear analysis is that 
a„ and t, must be computed in the deformed local region configuration. Having computed the 
normal and tangential stresses along boundaries CD and EF (see Figure 23), the collocation 
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STRUCTURE 

SUBSTRUCTURE 



Figure 31. Transverse Deformation, v, for Structure and Substructure Models, Nonlinear Analysis 
at 10 psi Pressure. 




procedure outlined in section 2.4 has the same format for the two types of analyses. For the 
loading conditions considered in the study, boundary DE of the local region is subjected to 
applied pressure, p. As stated earlier, this pressure is negligible, in the collocation procedure 
context, as compared to the stresses generated by the pressure (i.e., on boundaries CD and 
EF). 

Other restrictions in the present analysis which are a result of the assumptions used in the 
the elasticity solution relate to the uniformity of response in the local region. These re- 
strictions are: 

(a) The stresses do not vary with z. 

(b) e 0 , k, , and k,,, are uniform throughout the local region. 

For the skin-stiffener geometry, the loading conditions, and the material constitutive law used, 
these conditions are met to within a 10 % (or less) variation. With a 10 % variation, the re- 
strictions imposed by the elasticity solution are not perfectly satisfied. However, they are suf- 
ficiently close from an engineering view point. 

To utilized the convergence information established in section 2.6, the stress calculations in 
the local region using EAL solid elements are compared to the stress calculations using the 
PE2D finite-element. It was hoped that on the boundary of the local region, the EAL solid ele- 
ments would yield the same stresses as the PE2D elements. If the data compared favorably, 
then the convergence studies of Chapter 2, using PE2D, could be used to imply convergence 
of the EAL analysis. To facilitate this comparison, a case which was studied previously using 
PE2D was used. This is the case of bending illustrated in Figure 4. The results of one such 

comparison is shown in Figure 32. The comparison is for the coarse mesh Figure 4.a. As be- 

Mt 

fore, the interface stresses are normalized by Both the stresses in the skin and flange, 
at the interface, as computed by EAL and PE2D are shown. Further, the elasticity solution as 
derived from the collocation of stress data generated by EAL and generated by PE2D analyses 
are also plotted. Apparent from the figure is that the two approaches produce identical 
elasticity solutions. On the other hand, there are some differences with regard to the fmite- 
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Figure 32. Skin-Stiffener Stresses, a x , a , and t in the Local Region; EAL Verses PE2D. 



element computed stresses within the localized region. This is particularly apparent for . 
However, this is not surprising, since the EAL solid elements have constant shear stress 
through the element cross-section. Hence such an element will have a more difficulty con- 
forming to the large stress gradients in the localized region. Nevertheless, PE2D and EAL 
produce almost identical stresses along the local region boundaries atx/t, = ± 3. This in turn 
produces the perfect agreement between the two elasticity solutions. Based on the data pre- 
sented in section 2.6 and comparisons like the one shown in Figure 32, it was concluded that 
the use of EAL leads to an accurate local elasticity analysis. 

In the last part of this chapter the local analysis of the skin-stiffener interface in the pressure 
loaded plate is presented. This presentation completes the analysis procedure depicted in 
Figure 24, starting with the structural analysis of Figure 27, proceeding to the substructural 
analysis of Figure 30, and ending with results similar to Figure 32, but for a pressure loaded 
plate. In the following figures, the stresses are normalized by the applied pressure, p, and 
distance have been normalized by the flange half-width, b, or by the skin thickness, t,. 

The stresses along the entire skin-flange interface are shown in Figure 33. Figure 34 fo- 
cuses on the distribution of these same stress components in the local flange termination re- 
gion. It should be kept in mind that the loading is actually a 10 psi vacuum applied to the side 
of the plate without the stiffener. Thus the plate out-of-plane deformation is in the -v direction 
(see Figure 24). The flange termination angle is 90°. Figure 33 indicates that along the entire 
interface, on either side, the stress a x is positive. This despite the fact that the interface is 
on the top side of the skin and, in the sense of a linear analysis, subjected to compressive 
bending stresses. The positive a x is attributed to the tensile membrane force that develops 
as a result of the large out-of-plane deformations of the plate. 

The peeling, <r y , is nonzero at the flange termination point, the local region, and under the 
stiffener web. Since the stiffener web has much more resistance to out-of-plane deformation, 
the skin has a tendency to pull away at this locality as well as in the flange termination region. 
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The shear stress is nonzero along the entire skin-flange interface, but it is constant in 
magnitude away from the flange termination point, and away from the web. 

This concludes the discussion as related to the analytical developments and their applica- 
tion to actual stiffened-skin composite structures. In the next chapter the structure- 
substructure-local analysis procedure will be used to investigate the influence of geometric 
nonlinearities on skin-stiffener interface stresses. Following that is study of the influence of 
stiffener parameters on the state of interface 
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4. Results 


With the methodology for computing skin-stiffener interface stresses verified as regards 
convergence and fidelity, the influence of various stiffener parameters on the skin-stiffener 
interface stresses can be evaluated. In addition, the importance of geometric nonlinearities 
can be evaluated, and an assessment of the error incurred by not including these effects can 
be made. To that end, this chapter begins by evaluating the influence of geometric nonline- 
arities on the interface stresses. This is followed by the determination of the effects of stiffener 
parameters on the state of these stresses. In the last part of this chapter, criterion by which 
various stiffener designs may be evaluated for their tendency to separate, using the current 
analysis procedure, are presented. The results presented in the following sections are for a 
plate clamped on all four edges with a single stiffener, depicted in Figure 22. The plates are 
subjected to three levels of transverses pressure, 1, 10, and 20 psi, respectively. The 10 psi 
pressure level represents the operating aircraft fuselage pressure, while the 20 psi pressure 
represents the design level. The 1 psi pressure level is used since the linear and nonlinear 
analyses would produce similar results at this pressure level, and linear results for the other 
pressure levels can be scaled from the 1 psi analysis. 



4.1 The influence of Geometric Nonlinearity 


To a great extent, the importance of incorporating geometric nonlinearities in the skin- 
stiffener interface stress analysis will depend on how flexible the structure is. For stiffened 
aircraft structures which exhibit out-of-plane deformations on the order of magnitude of 2 to 
4 times the skin thickness, geometrically nonlinear effects can be significant. As was evident 
from Figure 26, the skin-stiffener configuration being studied did in fact experience such de- 
formation levels. 

To study the influence of geometric nonlinearities on the interface stresses, a baseline 
skin-stiffener configuration is considered. The three pressure levels are used, and the linear 
analysis case is also considered for comparison. The baseline stiffener design has a web 
height of h w = 1.5 in., a flange width of 2b = 1.5 in., and a flange thickness oft, = 0.04 in. 
For this configuration the flange and the web are constructed of a quasi-isotropic laminate 
with ( ± 45/0/90), layup sequence. The skin is orthotropic with a laminate layup of 
( ± 45/90;,),, and skin thickness oft, = 0.04 in. 

Figure 35 illustrates the skin-stiffener interface stresses, o*, o y , and along the entire 
flange-skin interface length. Each portion of Figure 35 illustrates four relations, three of which 
represent nonlinear analyses for 1, 10, and 20 psi pressures, and one that corresponds to the 
linear analysis. Each relation is determined using a combination of stress data generated by 
an elasticity analysis in the local region, and an EAL finite-element analysis outside of this 
locality. In Figure 35 and all other figures which display stress data, stresses are normalized 
by the applied pressure. With this normalization the linear analysis is independent of the ap- 
plied pressure level. The distance, x, from the flange termination vertex is normalized by b, 
the flange half width, or by t„ the skin thickness. Clearly evident from the figure is the non- 
linear interaction between the level of applied pressure and the magnitude of the skin-stiffener 
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interface stresses. Focusing on the peel stress, for example, the linear and the 1 psi nonlinear 
analyses produce stresses of similar magnitude. However, the nonlinear analysis for 10 psi 
and the linear analysis show significantly different results. The 20 psi nonlinear analysis de- 
viates even more from the linear analysis. In all cases and for all stress components, the large 
stress gradients appears to be confined primarily to the local flange termination region. Under 
the web there is a local maximum in the peel stress and the maximum value tends to de- 
crease slightly with increases in pressure. Relations like the one shown in Figure 35 indicate 
that load transfer from the skin to the stiffener is confined to either a very short portion of the 
flange width (something on the order of b/4), or directly under the web. In Figure 36 attention 
is shifted to the interface stress distribution in the flange termination area. As in the previous 
figure, there are four relations in each portion of the figure, one for the linear analysis and the 
other three for nonlinear analyses. Only the interface peeling, o y , and shearing, t^, stresses 
are shown in this figure. The most significant point to emerge from Figures 35 and 36 is a clear 
illustration of the nonlinear interaction between the applied pressure and skin-stiffener inter- 
face stresses. If a linear analysis is used, the normalized interface stress variation throughout 
the flange termination region would be independent of the applied pressure. However, with 
geometric nonlinearities included, significant changes in the characteristics of the stress dis- 
tributions are observed. In a geometrically nonlinear analysis an increase in pressure tends 
to flatten the peeling stress relation. Increases in pressure lead to a shift in the high stress 
gradients in both the peeling and shearing stress relations toward the flange termination 
point. As a result, it appears that increases in pressure lead to a reduction in the area under 
the two stress relations. The possible significance of such behavior in the stiffener separation 
context will be discussed latter. 

As noted above, the variation in o y throughout the local interface region is markedly dif- 
ferent between the 1 and 20 psi pressure levels. This may be understood by considering Fig- 
ure 37. In this figure the stress component a, along the collocation boundaries CD and EF (see 
Figure 23) is plotted as a function of the distance through the thickness, y. Boundary EF cuts 
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through the skin at x = - 3t„ whereas, boundary CD cuts through the skin and the flange at 
x = + 3t,. it should be noted that on these boundaries a x = cr n and = -r,, and that both ct„ 
and t, are collocated. Nevertheless, a n is about two orders of magnitude larger than x t , and 
therefore has the most significant influence on the local elasticity solution. As in previous fig- 
ures there are four relations on each plot, one for the linear analysis and three for the non- 
linear analyses at the 1, 10, 20 psi pressure levels. Concentrating for the moment on the 
stress in the skin, Figure 37.a, the most important point to emerge is that at a low applied 
pressure, the skin is loaded primarily in bending. As a matter of fact, in a linear analysis of 
the stiffened plate which is loaded by pressure, the skin should be loaded in pure bending. 
This is indeed evident from the line depicting the linear analysis results. For a linear analysis 
at higher applied plate pressures, the relation would remain one of pure bending. However, 
due to the geometric nonlinearities, the stresses in the skin change from a state which is 
nearly pure bending, at low pressure, to a state at the high pressure level which is dominated 
by stretching. At 20 psi the skin is entirely in tension. Turning to the skin-flange combination, 
the normal stress distribution through the skin and the flange at x = + 3t, (Figure 37.b) shows 
more complex characteristics. This is due to the abrupt change in material properties in going 
from the skin to the flange. Nevertheless, <y K becomes purely tensile at the higher pressure 
level. Since the local region must remain in force equilibrium, this is expected, the stress 
distribution in the skin-flange region balancing out the stresses in the skin away from the 
flange. The above data suggest that the observed differences in the peeling stress distribution 
for the 1 and 20 psi pressure levels are related to the state of inplane loading. That is, at low 
pressure levels, for which the plate is loaded primarily in bending, the peel stress distribution 
is high at the flange termination and reverses to become negative away from the terminus. 
At high pressure levels, for which the plate is loaded primarily in inplane stretching, the peel 
stress is reduced at the flange terminus and as a result does not experience as large a re- 
versal away from the terminus, i.e., the peel stress distribution flattens with distance from the 
terminus. 
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To facilitate a better demonstration of this important nonlinear interaction between skin- 
stiffener interface stresses and applied pressure, stress eigenfactors (SEF) are used. These 
stress eigenfactors are defined as: 

lim 

K v = x ^ a v (x, 0; L,) , [70.a] 

x - 0 

lim 

K xy = x ^ t„(x, 0; , [70.b] 

x - 0 

where K y and are the peeling and shearing stress eigenfactors, respectively, and X, is the 
first eigenvalue. As noted previously, values of - 1 £ Re(X.,) £ 0 will lead to singular stress 
characteristics near the flange termination vertex. It should be noted that in the field of frac- 
ture mechanics the stress eigenfactor is commonly known as the stress intensity factor. The 
stress eigenfactor is a measure of how rapidly the stresses become unbounded as the flange 
termination point is approached. The factor can be used as a measure of the severity of the 
interaction between the skin and the stiffener. Figure 38 shows the skin-stiffener interface 
peeling and shearing stress eigenfactors as a function of the applied plate pressure. The 
stress eigenfactors have been normalized by the peeling stress eigenfactor for a 1 psi linear 
analysis. The influence of nonlinear effects is clearly illustrated in the figure. For a linear 
analysis the stress eigenfactor would simply be in proportion to the applied pressure. How- 
ever, it is clear that inclusion of geometric nonlinearities results in a different behavior. Spe- 
cifically, the stress eigenfactors increase slower than the applied pressure. This is consistent 
with the finding of Figure 36 and both figures point to the need for considering geometric 
nonlinearities in the study of skin-stiffener interaction. More importantly, it would appear that 
when conducting failure analyses of such structures, it should be recognized that a doubling 
of applied pressure does not result in a doubling of interface stresses and to assume so would 
be in error. 
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Finally, the influence of geometric nonlinearities on the out-of-plane deformations of the 
baseline plate configuration are considered in Figure 39. The figure consists of two portions: 
one depicting the out-of-plane deformation along the symmetry line (x,0,0) and the other por- 
tion the deformation under the web, colinear to the stiffener, along the symmetry line (0,0, z). 
Each portion of the figure encompasses four relations, one representing a linear analysis at 
1 psi pressure, and the other three corresponding to the nonlinear analyses at 1, 10, and 20 
psi applied pressure. The out-of-plane displacement, v, is normalized by the skin thickness, 
t,. The distance from the plate center (in the directions x or z) is normalized by the flange 
half-width, b. Visible from Figure 39 is the large pillowing effect, the skin away from the flange 
experiencing much larger out-of-plane deformations than the flange. This pillowing effect be- 
comes more important with increases in pressure. As a result, geometric noniinearities have 
an important role in determining the character of the pillowing. Specifically, for a linear 
analysis at 10 psi pressure, the maximum skin deflection would be about 7t,, while the stiffener 
deflection would be about 2t,. The nonlinear analysis at the same pressure level indicates that 
the maximum skin deflection is only 2.4t„ while the stiffener deflection is 1.4t,. Therefore it is 
clear that geometrically nonlinear effects actually reduce pillowing. The pillowing effect shown 
here was observed experimentally in [10], 

This concludes the discussion of the importance of geometric nonlinearities in the analysis 
of skin-stiffener interface stresses. It has been shown that neglecting these effects may lead 
to large errors in the calculation of plate response, particularly interface stresses. Next, at- 
tention is given to an examination of the influence of stiffener geometry and material proper- 
ties on skin-stiffener interface stresses. 


4.2.1 Stiffener Parametric Study Results 






To establish the usefulness of the analysis method in the design of stiffened composite 
structures, eight different stiffener configurations are studied. Table 2 lists the various stiffener 
designs considered, and their designation, geometry, and lamination sequence. In the fol- 
lowing sections each stiffener configuration is referred to by a nickname (second column in 
Table 2). This nickname signifies what is unique about this stiffener relative to the baseline 
configuration. The baseline configuration, configuration A, is always referred to by the nick- 
name 'baseline*. The baseline stiffener design was defined earlier but for the reference, it has 
a web height h w = 1.5 in., flange width 2b = 1.5 in., and flange thickness t, = 0.04 in. For 
baseline, the flange and the web are constructed of a quasi-isotropic laminate with 
(±45/0/90), layup sequence. The skin is orthotropic with a laminate layup of (±45/90 2 ) t . In 
comparison to the baseline, the nickname "soft flange' was chosen for configuration B since 
the flange layup resulted in transverse flange modulus E x = 3.83 Msi relative to E„ = 8.05 
Msi of the baseline configuration. All other stiffener nickname designations follow the same 
logic. Lamina and laminate material properties are given in Appendix D. The structural anal- 
ysis is conducted using the mesh shown in Figure 27. It should be noted that tapering of the 
flange will lead to a slight reduction in the moment of inertia of the flange relative to the case 
of a 90° flange termination angle. However, the structural analysis model is insensitive to the 
influence of the reduced moment of inertia. For this reason configurations A, B, D, and F are 
analysed using the finite-element discretization model for the baseline flange. Similarly, con- 
figurations E, G, H, and I are analysed using the finite-element discretization model for the 
thick flange. However, the substructural finite-element model for 90° and 15° flange termi- 
nation angles are not the same, the difference reflecting the flange angle. The difference is 
similar to the difference between Figures 4 and 5. Based on the results which were presented 
in Figures 18 and 19, it is possible to use only one substructural finite-element analysis, for 
example, the one for the 90° flange angle, to produce two local elasticity solutions, one for 
90° flange termination angle and another one for 15° flange termination angle. However, 15° 
is considered a shallow angle and it is not clear that one substructural analysis is all that is 
necessary. Since an important part of the current investigation is related to the method de- 



Table 2. Stiffener Configurations. 



Thick, Soft, & Tapered Flange 1.50 0.08 15° ( ± 45/0 4 / ± 45) 


























































velopmenl and verification, two separate substructure analyses are used. In a latter part of 
this section data is presented to support the assertion that a single analysis would lead to 
acceptably accurate results, but not as accurate as for less shallow angles, say 45°. 

The results in the following section are presented in a form of skin-stiffener interface peeling 
and shearing stresses throughout the localized region. Each figure represents a comparison 
of the stresses between the baseline stiffener configuration and one of the other eight stiffener 
designs. Stress data are given for 1 and 20 psi applied pressure and are based on a nonlinear 
analysis. The low pressure represents the level at which linear and nonlinear analyses give 
similar results. The high pressure level represents a pressure where geometrically nonlinear 
effects become significant. The figures are useful in demonstrating qualitatively how interface 
stresses are affected by the various stiffener parameters. However, they may be less helpful 
quantitatively since they provide stress distributions rather than one numerical value which 
can be applied to design, failure, or optimization procedures. The calculation of a numerical 
value is done in the last part of this chapter, where the performance of each stiffener config- 
uration is evaluated against the baseline stiffener design. 

Figure 40 shows a comparison of skin-stiffener interface peeling and shearing stresses 
between the baseline and the "soft flange" configurations. It appears from the figure that re- 
ducing the transverse flange modulus, E x , more than 50% has little effect on the peeling 
stress. It does tend to slightly reduce the shearing stress near the flange vertex. Figure 41 
details the interface peeling and shearing stresses for both the baseline and the "short web" 
configurations. Evident from the figure is that shortening the web leads to a significant re- 
duction in both the peeling and shearing stresses. At the high pressure, the peeling stress is 
essentially zero, while the shearing stress, although reduced, is still measurable. Since the 
softer stiffener configuration resulting from the short web leads to lower bending gradients in 
the flange termination region, the reduction in the peeling stress is somewhat expected. 











Figure 42 illustrates the character of the peeling and shearing stresses for the baseline and 
"tapered flange" stiffener configurations. Tapering the flange seems to have a considerable 
effect on the interface peeling stress distribution. At both low and high pressures it appears 
that tapering tends to distribute the stress as a tensile stress over a larger portion of the 
skin-stiffener interface, rather than having a stress reversal accompanied by a steep gradient. 
This, in turn, may reduce the separation tendency of this particular stiffener design. Tapering 
of the flange seems to have a smaller influence on the interface shearing stress distribution. 
In general it tends to slightly increase for both low and high pressure levels. Figure 43 
demonstrates the influence of thickening the flange on the skin-stiffener peeling and shearing 
interface stresses. Generally speaking, thickening the flange tends to increase the interface 
peeling stress, increasing the gradient near the flange terminus and increasing the reversal 
away from this point. On the other hand, thickening the flange has a mixed effect on shearing 
stress. At high pressure the shear stress increases, while at low pressure the shear stress 
increases or decreases, depending on spatial location. 

Figure 44 delineates the effect of softening and tapering the flange at the same time. An 
examination of Figure 42 and Figure 44 reveals that the peeling stress distributions for the 
tapered flange configuration and for the soft and tapered flange configuration are very similar. 
However, in the latter case the stresses appear to be lower, reflecting the influence of the 
combined changes to the baseline case. The shearing stress distribution seems to be lower, 
as compared to the baseline design, near the vertex and slightly higher away from the vertex. 

Figure 45 displays a comparison between the baseline configuration and one for which the 
flange was thickened and tapered. If the results of Figure 45 are evaluated in the context of 
the results of Figure 43, it may be concluded that tapering the thick flange improves the dis- 
tribution of the peeling stress along the interface, but has a smaller effect on the shearing 


stress distribution. 




















Figure 46 depicts the influence of both thickening and softening the flange on the peeling 
and shearing interface stresses. Again, if one examines the results of Figure 46 in light of the 
data which were presented in Figure 43, it appears that softening the thick flange has little 
influence on the state of interface stresses as compared to the "thick flange" design. There- 
fore, the discussion offered for Figure 43 is also relevant for Figure 46. These results are also 
consistent with the results of Figure 40, where the baseline stiffener configuration was evalu- 
ated relative to a "soft flange" stiffener design of the same geometry. Finally, Figure 47 dem- 
onstrates the effect of thickening, softening, and tapering the flange on skin-stiffener interface 
peeling and shearing stresses. The following observations are made in connection with 
Figure 47: Softening the tapered thick flange does little to change the peeling stress distrib- 
ution in the local region (compare Figure 45 with Figure 47). However, tapering the soft thick 
flange does substantially reduce the peel stress (compare Figure 46 with Figure 47). Finally, 
tapering the soft thick flange substantially reduces x^. Based on the results of Figure 47 and 
the other figures, it can be stated that the single most important influence on the interface 
stresses is flange tapering. Softening the flange and reducing the thickness have considerably 
less influence. 

At this point it is of interest to determined if the use of a single substructural analysis is 
sufficient for obtaining a local elasticity solution for both 90° and 15° flange termination angles. 
In section 2.6 it was shown that for the 90° and 45° flange termination angles, the local 
elasticity solution could be obtained from one substructural (or global) finite-element analysis 
for either 90° or 45° flange termination angles. However, some concern was raised as to 
whether these results are applicable for a very shallow flange termination angle. The results 
illustrated in Figure 48 provide insight into this. Figure 48 replicates the previous study of the 
45° and 90° flange termination angles (Figures 18 and 19) but considers 15° and 90° angles 
instead. In Figure 48 the the pressure loaded stiffened plate and clamped boundary conditions 
used in the above parameter study are considered. The two stress components of interest, 
Gy and T,y, are illustrated. The solid line in each portion of the figure represents the local 
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Figure 47. Influence of Thickening, Softening, and Tapering the Flange on the Skin-Stiffener 
Interface Peeling and Shearing Stresses. 










elasticity solution obtained from the collocation of stress data from the 90° flange termination 
angle substructure finite-element analysis. On each portion of the figure there are two solid 
line relations; one illustrates the local elasticity solution for a 90° flange, and the other illus- 
trates the solution for a 15° flange. These curves are labeled "90° flange" and "15° flange", 
respectively. Similarly, the dashed lines on each portion of the figure represent the local 
elasticity solution obtained from the collocation of stress data from the 15° flange termination 
angle substructure finite-element analysis. There are two dashed line relations on each por- 
tion of the figure; one representing the local elasticity solution for the 90° flange, the other 
representing the local elasticity solution for the 15° flange. If the two curves which are labeled 
"90° flange" (i.e., the dashed and the solid lines) coincide, it is possible to use either the 90° 
or 15° finite-element analysis results to produce the same 90° flange elasticity solution. Simi- 
larly, if the two curves labeled "15° flange" coincide, the use of either 90° or 15° flange ge- 
ometry finite-element analysis results to produce a 15° flange elasticity solution is justified. In 
general, the results of Figure 48 indicate that the level of agreement between the dashed and 
the solid Ijnes for the same local elasticity solution produced by two different substructural 
analyses is not as good as the one observed in Figures 18 and 19. From a stiffener failure 
analysis point of view (discussed in the next section) the results of Figure 48 may be close 
enough. The worsening in agreement between the two approaches is directly related to the 
shallow flange termination angle, 15°. It is reasonable to expect that at very shallow flange 
termination angles, modeling the flange tapered portion in the substructural analysis is im- 
portant. In addition, the results of Figures 18 and 19 are based on a linear analysis, whereas 
the present results are based on a geometrically nonlinear analysis. It is possible that the 
interaction between geometry and flange-skin interface response are more significant in the 
nonlinear range. 

Lastly, a significant point of interest is the out-of-plane deformations observed in the various 
stiffened plate configurations. Figure 49 illustrates the plate deformation v at the center (z = 0) 
as a function of transverse distance, x, for four stiffener configurations. These configurations 






are: the baseline case, the short web case, the thick flange case, and the soft flange case. 
The figure demonstrates some important aspects of skin-stiffener interaction. Initially, at low 
pressure levels, all skin-stiffener configurations show a pillowing effect. At high pressure level, 
the short web stiffener configuration does not pillow. The other 3 stiffener configurations, 
however, exhibit pillowing at high pressure levels. However, there are differences among 
these three cases. The thicker flange stiffener geometry results in nearly a rigid body trans- 
lation downward at the flange, i.e., v in the flange is practically independent of x. The soft 
flange, on the other hand, bends, its downward deflection increasing with increasing x. The 
baseline case represents a deformation characteristic between these two extremes. It should 
also be noted that the plate with the soft flange deflects less than the plate with the baseline 
flange. This is a consequence of the fact that in the soft flange configuration the 90° fibers in 
the baseline flange are converted into 0° fibers. Since it is assumed that the flange and web 
are of the same laminate, the switching of the 90° fibers to 0° fibers stiffens the web, thereby 
decreasing the overall deflection of the plate. 


4.2.2 Stiffener Design Performance Evaluation 


In the above, the influence of stiffener design parameters on skin stiffener interaction have 
been examined from a qualitative point of view. In the following section, these effects are 
studied quantitatively. Interest here is in presenting numerical values by which a particular 
skin-stiffener design may be judged in terms of its tendency to experience skin-stiffener sep- 
aration failure. The results are presented relative to the baseline configuration. In the results 
to be presented, two stress parameters are used, the stress eigenfactor (SEF), mentioned 
previously, and the average stress factor (ASF). 



The peeling, Ky, and shearing, stress eigenfactors, were defined by eq. 70.a and 70.b, 
respectively. In the present investigation the SEF concept is used as a measure of the severity 
in interaction between the stiffener and the skin, rather than as a failure criteria. A stiffener 
design which leads to a high SEF, relative to a baseline configuration, will be judged as having 
more tendency to separate. Similarly, a stiffener configuration which produces a lower SEF, 
relative to a baseline design, will be judged as having less tendency to separate. For example, 
if tapering the flange of the baseline configuration leads to lowering the SEF, then reducing the 
stiffener termination angle will be judged as a design improvement for the stiffener. Here care 
must be taken in evaluating the results. Changing one stiffener parameter my lead to a re- 
duction in Ky but an increase in K,y, or vice versa. This raises the question of which SEF is 
more significant in the initiation of skin stiffener separation K y or K^. Here it can only be said 
that epoxy resins and adhesives usually exhibit higher shearing strength relative to peeling 
strength. Therefore, it is reasonable to speculate that K y will be more important in determin- 
ing the tendency for skin stiffener separation. 


While the concept of a stress eigenfactor has been used successfully in the field of fracture 
mechanics (i.e., stress intensity factor) in isotropic materials as a parameter by which the in- 
itiation of crack growth and fracture can be determined, this concept has proven to be less 
useful in composite materials. This may be related to the heterogeneity of the material and the 
complex nature of crack propagation. In the skin-stiffener separation initiation context, there 
is no evidence for the existence of crack at the flange terminus prior to failure. Therefore, the 
use of the average stress factor in determining skin-stiffener interface strength may be more 
appropriate. For the present problem, the average stress factors are defined by; 


S„ = 
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where is the i'th eigenvalue in the truncated eigenfunction expansion. The parameter a 0 is 
the distance of integration along the flange interface from the flange terminus. Failure is based 
on critical values for S" and S" for a particular material. These must be determined exper- 
imentally. The average stress criteria was first proposed for use in composite structure ap- 
plications in [45], In that study the criteria was applied to the calculation of laminates strength 
which contained holes and notches. More recently this concept was used in the investigation 
of the onset of delamination in composite laminates [46], In that study the average stress cri- 
teria produced fairly good agreement with experimental results. The investigators postulated 
that the onset of delamination was strictly controlled by interlaminar transverse strength and 

Q 

therefore used S y only. The integration length was taken as -p = 1 (t, being the lamina 
thickness) and S" was equal to the transverse strength of the composite. Due to the lack of 
experimental data, the average stress criteria in the present investigation is used as an ad- 
ditional parameter by which the different stiffener designs are evaluated. It is anticipated that 
both the SEF and the ASF will produce qualitatively similar results. 

Finally, it is recognized that skin-stiffener designs cannot be judged solely by their tendency 
to separate. Therefore, in addition to the stress parameters, the maximum skin and the 
center-of-plate out-of-plane deformation are given. This is to point out that a particular 
stiffener design may lead to a 50% reduction in K, and but an increase of 2 to 3 times in 

the out-of-plane deflection. Such designs would be desirable from the skin-stiffener separation 
stand point, but may be unacceptable from the structural deformation point of view. 

Table 3 presents quantitative data relevant to the stiffener parametric study. The data given 
in the table is based on a 1 and a 20 psi pressure nonlinear analysis. The first column lists the 
stiffener configurations by their nicknames. The second column gives the two pressure levels. 
Columns 3 and 4 give the peeling and shearing stress eigenfactors, K y and K^. Columns 5 and 
6 provide the peeling and shearing average stress factors, S y and S^,, for an integration dis- 
tance of aj\, = 0.0625. This value of a 0 corresponds to one lamina thickness for the materials 
considered here. Finally, columns 7 and 8 list values for the maximum skin and the center- 



> 


M 

0) 

AC 

>* 

-o 

3 

(/) 

o 

u 

<D 

E 


£ 

L. 

4 ) 

C 

S 

‘-5 

W 


CO 

0) 

43 

|2 


>“ 

CD in 
CO 04 

h- co 

CM CO 
CO CM 

xt xt 

II 

CD x- 
xt 
m co 

cd m 

oo CM 
h- oo 

CO X- 

r^- xr 
m cd 

05 XT 
XT x- 

xr in 

05 XT 
XT x- 

xr in 


o o 

in cm 


oo 

o o 

o o 

o o 

d o 

M 

> 

T- Tf 
x- m 

o 05 

CD Tf 
CM 05 
t*- xt 

II 

CD O 
05 CO 
00 00 

X- Tf 

x- in 

O 05 

CD o 

05 CO 

00 oo 

CO 1^ 
05 05 

oo h- 

co h- 

05 05 
00 h- 


x- o 

T~ T* 


o o 

x- o 

o o 

o o 

o o 

fr 

CO 

00 N- 
xT TT 
h- h- 

xt o 

05 05 
in m 

05 CM 
CO CD 
CD CD 

CD CO 
xr- CM 
x- TT 

CO h- 
xf CD 
CO CO 

CM 05 
00 I s - 
XT CD 

h- 00 
05 CD 
oo x- 

cd 
in r^. 
CM co 

o o 

o o 

o o 

T~ X** 

o o 

o o 

O x- 

o o 

>» 

CO 

O) TT 
CM CM 
05 05 

CO Tf 

co xt 
in cm 

CM O 
r^. x- 
-r- CM 

xr* CO 
CM 00 
CM CD 

Q CD 
O CM 

oo CD 
CM x- 
x- CM 

05 05 

in oo 
cm 

XT 05 
h- CM 
O x- 


o o 

o o 

o o 

T- T— 

o d 

o o 

x— x— 

o o 

fr 

* 

O CD 
CM x- 
CO oo 

co o 

t- 

in m 

in co 
in 

O T- 

x- CM 
CD CM 
t- m 

in 
oo r^ 
XT 05 

CM O 
x— CO 
x- CD 

in xr 

CO CO 
O CO 

XT CO 
CO oo 
CD O 


o o 

o o 

xt xt 

r* x— 

CO CO 

CO XT 

x— x— 

CM xr 

* 

I s - CO 

S3 

co o 

i**- x— 

in in 

co 
o r*- 

CO CO 

x- CM 
CD CM 

x-; in 

X- oo 
X- CM 

in cd 

XT O 
CO O) 
XT to 

CM 05 
CM CD 

oo 

cd in 

X- XT 
XT CD 


T* T* 

d d 

o d 

r- x- 

O O 

d O 

x- X— 

o o 

ressure 

(psi) 

1 

20 

^ 04 

i 

20 

1 

20 

i 

20 

1 

20 

1 

20 

_ o 

^ CM 

a. 









Stiffener 

Configuration 

0) 

O) 

c 

CO 

lZ 

% 

CO 

n 

0) 

5 

t: 

o 

JO 

CO 

Tapered Flange 

0) 

o> 

c 

« 

EZ 

JX 

G 

is 

Soft & Tapered 
Flange 

Thick & Tapered 
Flange 

Thick & Soft 
Flange 

Thick, Soft, 

& Tapered Flange 


CD 

d 

d 

ui 

lZ 

6 

x 

±J 
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of-plate out-of-plane deflections, v, and v f , respectively. All values which are given in Table 3 
are normalized by the values calculated for the baseline configuration. Table 3 presents an 
enormous amount of data that can be compared in ways that are too numerous to be done 
here. However, two types of comparisons can be done, namely comparison resulting from 
softening, tapering, and thickening the baseline flange, and softening and tapering the thick 
flange. These comparisons will be made in the following paragraphs. 

Examination of the SEF data lead to the following observations: softening the flange (i.e., 
reducing the transverse modulus EJ tends to increase K y , in both the baseline and thick 
flanges. This increase is more significant in the thick flange, as can be seen by a comparison 
between the SEF for the "thick flange" and the "thick and soft flange" configurations. On the 
other hand, softening the flange leads to a reduction in in both the baseline and thick 
flanges. For the soft flange, the reduction is 20% relative to the baseline case. Shortening 
the web, to half the height of the baseline configuration, results in drastic a decrease in both 
the peeling and shearing stress eigenfactors. Though not as dramatic, tapering the flange 
leads to a 30% or more reduction in Ky in the baseline case. For the thick flange configuration 
tapering leads to even greater reduction in Ky. However, this same design change leads to 
large increases in K v Thickening the flange brings about an increase of up to 50% in both K y 
and K«y at 20 psi pressure. This increase happens even though the values of the engineering 
constants are the same for the two flanges (see Appendix D). Softening the tapered flange 
helps reduce Ky and K^ by a moderate amount. This reduction amounts about 15% decrease 
in Ky and K,*, respectively, for the thick flange at 20 psi pressure. A smaller decrease is ob- 
served at a lower pressure and for the baseline flange configuration. On the other hand, ta- 
pering the soft flange substantially reduce Ky in both the baseline and thick flanges 
configurations. This design change brings about a 50% reduction in Ky, for the baseline flange 
at both 1 and 20 psi pressure. The same design change leads to about 65% reduction in Ky 
at 1 and 20 psi pressure in the thick flange configuration. However, tapering the soft flange 
leads to a drastic increase in the shearing stress eigenfactor in both the baseline and thick 


flange configurations. This increase is similar in magnitude to the increase due to tapering 
alone. Based on the SEF data presented, the following general comments can be made: (a) 
Softening the flange leads to slight increases in K y and moderate decreases in K^; (b) Taper- 
ing the flange brings about a drastic reduction in K y but a substantial increase in K v ; (c) 
Thickening the flange leads to up to a 50% increase in K y and K^; and (d) Shortening the web 
heights to half the baseline web height results in nearly 50% reduction in both K y and K v . 

The average stress factor produces results which, for the most part, are consistent with the 
SEF data. However, there are some inconsistencies.These are discussed next. The ASF, unlike 
the SEF, indicates that softening the flange leads to slight decrease in peel tendency. How- 
ever, the difference between the two measures is quite small, 5% increase in K y verses 8% 
decrease in S y . Another inconsistency relates to the values of and S^ for the "thick & soft 
flange" configuration at 1 psi pressure. For this configuration K ^ is 4% higher than the 
baseline design, whereas S^ is 10% lower. However, it should be noted that softening the 
thick flange produces the same trend in both the SEF and the ASF. That is, both K y and S y in- 
crease due to softening, whereas K v and S v decrease due to softening. Finally, in the SEF 
calculation, tapering the flange leads to an increase of 3 to 4 times in relative to the 
baseline configuration. The ASF shows a decrease in S^ relative to the baseline stiffener de- 
sign. As a matter of fact, all the stiffener configurations except the "thick flange* and the "thick 
& soft flange' designs produce S^ which is smaller than the baseline case. These results are 
consistent with the observations made in Figure 40 through Figure 47. The average stress, 
of course, will depend on the distance of integration, a 0 . For example, if the "tapered flange" 
configuration, Figure 42, is considered, it is obvious that extending the integration distance to, 
say, aj\ t = 1 would result in different values for both S y and S^. This is what makes the SEF 
so useful in evaluating the effects of various stiffener parameters on skin stiffener interaction. 
The SEF is independent of any objective parameter such as a 0 . On the other hand, the ASF 
may be useful in determining the skin stiffener separation failure initiation. This can only be 
achieved with experimentally determined value for S" and S" 



Finally, an examination of the out-of-plane deflections for the various stiffener designs 
shows that, relative to the baseline case, all configurations, except for the "short web’ design, 
give the same or smaller out-of-plane deflections for the center-of-plate and the skin. Although 
the 'short web' stiffener design leads to a substantial reduction in both the peeling and 
shearing stresses in the local region, it results in a large increase in the out-of-plane de- 
flections. 

In summary, if the stress data is considered together with the out-of-plane deflections, it 
may concluded that: (a) If a thick flange is used to increase structural rigidity, tapering and 
softening the flange will substantially reduce the risk of premature skin-stiffener separation; 
(b) Softening the flange (i.e., putting more 0° fibers in the z direction) leads to a decrease in 
the center-of-plate out-of-plane deflection but a relatively small gain as far as interface 
stresses are concerned; and (c) Tapering the flange results in a small or nonexistent effect 
on the out-of-plane deformation but substantially reduces the risk of premature skin-stiffener 
separation. 

This concludes the quantitative evaluation of the performance of stiffener design parame- 
ters. In the next chapter some conclusions are put forth and recommendations for future re- 
search, as related to skin-stiffener interaction, are discussed. 
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5. Conclusion and Recommendations 


5.7 Concluding remarks 


In the following chapter, the work conducted in the present study is summarized. In addition, 
some recommendations for future research are highlighted. 

As stated in the introduction, the objective of the investigation was to develop an analytical 
method by which skin-stiffener interface stresses could be accurately calculated. Particular 
attention was given to the flange termination region, a region where the stresses are high due 
to geometric and material discontinuity associated with this locality. Moreover, the stresses 
in this region are generally difficult to compute and so extra attention was given to this area 
of the interface. Furthermore, since stiffened panels are most commonly designed to operate 
at the postbuckling range, the analysis had to incorporate geometrically nonlinear effects. 
Considering the above, the following objectives were set: 

(a) The analytical model should accurately represent the state of stress near the 
point of geometric and material discontinuity. 
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(b) The analysis should include geometric nonlinearities. 

(c) The procedure should be applicable to general stiffened composite panels. 

(d) The model should be sensitive to various skin and stiffener design parameters, 
such as stiffener geometry, and stiffener/skin material architecture. 

The method developed addressed all of the above objectives. In Chapter 2 the first objective 
was addressed. This led to the development of a local-global analysis procedure, later called 
the local-substructure analysis procedure. The developments in this chapter were for ge- 
ometrically linear analyses. The local elasticity analysis was based on the eigenvalue expan- 
sion of the stress function. The eigenvalue expansion was applicable in the flange termination 
region and as a result, the stresses were known to within a set of arbitrary, but unknown, co- 
efficients which were associated with the eigenvalues. For a particular skin-stiffener problem, 
these coefficients were determined from the stresses obtained from a combination of a global 
finite-element analysis of the entire skin-stiffener cross-section and a collocation scheme. 
Once these constants were known, the stresses in the localized region could be uniquely de- 
termined. The local elasticity solution provided a rigorous solution which accurately charac- 
terized the material and geometric discontinuities associated with the flange termination 
region. Away from this region the global finite-element analysis produced accurate skin- 
stiffener interface stresses. The last part of Chapter 2 was devoted to the study of conver- 
gence and accuracy of the local elasticity solution. The accuracy issue was investigated by 
comparing the local elasticity solution with finite-element results for which the mesh was re- 
fined twice. The use of finite elements for comparison was due to the lack of another analysis 
of this particular problem. In general, it was found that the local elasticity solution scheme 
produced a very accurate interface stress representation in the flange termination region. 
Convergence of the local elasticity solution was studied by varying the number of eigenvalues 
in the truncated eigenvalue expansion, and by varying the number of the collocation points 
on the local region boundary. Generally, it was found that the use of 10 to 15 eigenvalues, in 
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the eigenvalue expansion, and 100 collocation points resulted in a converged local elasticity 
solution. 

In Chapter 3 the second and third objectives were addressed. In the first part of the chapter 
the local elasticity solution was extended to include geometrically nonlinear effects. Here the 
condition of small strains but finite rotations was assumed. It was shown that under these 
conditions, in a local region which undergoes spatially uniform rotations, the local elasticity 
solution procedure developed in Chapter 2 was valid, provided it was applied in the deformed 
body configuration. This approach was applied to an actual stiffened composite plate structure 
where both geometrically linear and nonlinear analyses were considered. The analysis of the 
entire structure consisted of three steps: a structure analysis, a substructural analysis, and a 
local elasticity analysis. First a structural level analysis of the entire stiffened plate was con- 
ducted using finite-element plate elements. Next, a typical cross-section region in the stiffened 
plate structure was isolated and a three dimensional finite-element analysis of the cross- 
section was performed. Finally, the stresses from the substructural analysis were coupled with 
a boundary collocation scheme to produce a rigorous elasticity solution in the flange termi- 
nation region. 

Finally, in Chapter 4 both the third and the fourth objectives were addressed. In this chapter 
the influence of geometric nonlinearities on skin-stiffener interface stresses was evaluated. 
The interest here was focused on determining what error would be encountered if a geomet- 
rically linear, rather than a geometrically nonlinear analysis was used in the computation of 
skin-stiffener interface stresses. In general, it was found that in flexible stiffened skin struc- 
tures, which exhibit out-of-plane deformations on the order of magnitude of 2 to 4 times the 
skin thickness, geometrically nonlinear effects in the calculation of interface stresses are very 
important. That is, the use of geometrically linear analysis, rather then nonlinear analysis, can 
lead to considerable error in the computed interface stresses. In the last part of the chapter 
the influence of stiffener design parameters on the skin-stiffener interface stresses was 
studied. Both geometric and material stiffener parameters were considered. The stiffener 
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parametric study included eight different stiffener configurations. In general, it was found that 
the local elasticity analysis was sensitive to the various stiffener design changes and provided 
very useful information by which these designs could be evaluated against a baseline stiffener 
configuration. The findings suggest that tapering and softening the flange simultaneously 
tends to reduce interface peeling stress in the local flange termination region. However, these 
same design modifications tend to increase the shearing stress in this locality. Since it is 
believed that the stiffener separation failure is strongly related to the state of peeling stress 
at the flange termination region, the use of flange tapering and softening may prove benefici- 
ary in this respect. Finally, use was made of the peeling and shearing stress eigenfactors 
(SEF) in evaluating the performance of the various stiffener designs against a baseline con- 
figuration. Another stress parameter used was the average stress factor (ASF). Both the SIF 
and the ASF could be calculated in a closed form from the truncated eigenvalue expansion 
elasticity solution. 

As summarized above, all the objectives were addressed by the method developed in this 
study. The results demonstrated that this analytical tool is accurate and sensitive in solving 
skin-stiffener interface stress problems in actual stiffened composite skin aircraft structures. 


5.2 Recommendations for future research 


The method developed produced accurate and useful skin-stiffener interface stress calcu- 
lations for stiffened composite skin aircraft structures. However, the analytical model was 
based on certain assumptions which restricted its application to the most general skin- 
stiffener interface problems. In this section recommendations are made related to the capa- 
bility enhancement of the current model to a more general one. Other recommendations are 
related to analysis cost reduction, and to analytical and experimental correlation of stiffener 
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separation predictions. The recommendations are divided into two categories, one which re- 
lates to analytical development and the other is related to experimental considerations. 


5.2.1 Analytical Recommendations 

(a) In the present investigation the skin and stiffener were restricted to symmetric balanced 
laminates. This, in conjunction with the integrated laminate properties assumption, led to 
orthotropic material constitutive relations. As a result, material stretching and bending and/or 
twisting coupling effects were ignored. For most stiffened composite skin aircraft structures 
such a simplification is justified. However, there are some cases in which the use of unbal- 
anced or unsymmetric laminates for the skin and/or the stiffener is desirable. Hence, to render 
the present analytical tool more general, the local elasticity solution should be extended to 
include more general material constitutive law. It should be noted that for more general ma- 
terial characteristics, the two stress functions, F(x,y) and 'P(x.y) are coupled. However, in 
general, the solution approach is the same (see refs. 21 and 22). In that regard, it should be 
mentioned that for the cases studied here, and x„ were identically zero. 

(b) In the present investigation it was assumed that the stress components did not vary with 
z. This restriction is justified for most instances over short distance in the z direction where 
the local elasticity solution is applied. However, there are some cases in which the stresses, 
and in particular a z , do vary rapidly with z. In such cases the local elasticity solution can be 
modified to include a linear variation in a x with z (see ref. 22). The solution procedure to such 
a case is similar to the one presented here. 

(c) Finally, in the present analysis procedure use was made of the substructural analysis 
of the skin-stiffener cross-section. The primary function of this step was to generate the 
boundary conditions for the local elasticity analysis, as well as for the computation of skin- 


5. Conclusion and Recommendations 


135 



stiffener stresses outside of the local region. If the sole purpose of this step is only to generate 
boundary conditions, it would be advantageous from a cost standpoint to eliminate this step. 
One way by which the above may be achieved is by using laminated plate theory to compute 
the stresses on the local region boundaries using the inplane loading from the structural 
analysis. By doing so the local elasticity analysis can be incorporated as a post-processor in 
a structural analysis program. 


5.2.2 Experimental Recommendations 

(a) The lack of experimental data makes it impossible to evaluate the usefulness of the 
present analytical approach in determining the initiation of skin-stiffener separation failure. 
Therefore, it is recommended that stiffener-skin specimens configurations will be tested to 
failure under various loading conditions to obtain skin-stiffener interface strength, S‘ r and S" 

(b) Once critical values of the failure parameters are known, they can be used as input to 
make analytical failure predictions of actual skin-stiffener structures. For example, using the 
experimental devise describe in [10] for testing single stiffener pate configuration subjected 
to uniform transverse pressure, a comparison can be made between analytical prediction and 
experimental results. 


5. Conclusion and Recommendations 


136 


6. References 


1. Starnes, J.H., Jr. and Knight, N.F., Jr., ' Postbuckling behavior of selected flat stiffened 
graphite-epoxy panels loaded in compression," AIAA Paper No. 82-0777, presented at the 
AIAA/ASME/ASCE/AHS 23rd SDM Conference, 1982. 

2. Knight, N.F, Jr and Starnes, J.H., Jr., " Postbuckling behavior of selected curved stiffened 
graphite-epoxy panels loaded in compression," AIAA Paper No. 85-0768, presented at the 
AIAA/ASME/ASCE/AHS 26th SDM Conference, 1985. 

3. Knight, N.F., Jr. and Stroud, W.J., "Computational structural mechanics: A new activity at the 
NASA Langley Research Center," NASA TM 87612, Sept. 1985. 

4. Lotts, C.G., Green, W.H., McCleary, S.L., Knight, N.F., Jr., Paulson, S.S., and Gillian, R.E., ' 
Introduction to the Computational Structural Mechanics Testbed,", NASA TM 89096, Septem- 
ber 1987. 

5. Wang, J.T.S. and Biggers, S.B., ' Skin/stiffener interface stresses in composite stiffened 
panels," NASA Contract Report 172261, 1984. 

6. Goland, M. and Reissner, * The stresses in cemented joints," J. of App. Mech. , 1, 1 (1944) 
p. A-17. 

7. Yuceoglu, U. and Updike, D.P., " Stress analysis of bonded plates and joints," J. Eng. Mech. 
Div. , 106, (1980) p. 37. 


6. References 


137 



8. Williams, J.H., Jr., ' Stresses in adhesive between dissimilar adherents, ' J, Adhesion , 7, 
(1975) p. 97. 

9. Delale, F„ Erdogan, F., and Aydinoglu, M.N., ' Stresses in adhesively bonded joints: A 
closed-form solution,' J, of Composite Materials , 15, (1981) p. 249. 

10. Loup, D.C., Hyer, M.W. and Starnes, J.H., Jr., " Stiffened skin interaction in pressure loaded 
composite panels,' AIAA Paper No. 86-0917, presented at the AIAA/ASME/ASCE/AHS 27th 
SDM Conference. 

11. Williams, M.L, ' Stress singularity resulting from various boundary conditions in angular 
corner of plates in extension,' J. of App. Mech., 19, 4 (1952) p. 526. 

12. Gross, B., Srawley, J.E. and Brown, W.F., ' Stress intensity factors for single edge notch 
tension specimen by boundary collocation of stress function,' NASA TN D-2395 1964. 

13. Gross, B. and Srawley, J.E., ' Stress intensity factors for three points bend specimen by 
boundary collocation,' NASA TN D-3092 1965. 

14. Gross, B. and Srawley, J.E., ' Stress intensity factors by boundary collocation method for 
single-edge-notch specimens subjected to splitting forces," NASA TN D-3295 1966. 

15. Gross, B. and Srawley, J.E., ' Stress intensity factors for crackline-loaded edge-crack 
specimens,' NASA TN D-3820 1967. 

16. Gross, B. and Mendelson, A. ' Plane elastostatic analysis of v-notched plates,' Int. J. of 
Fracture Mech. , 8 (1972) p. 267. 

17. Carpenter, W.C., 'Calculation of fracture mechanics parameters for general corner," Int. 
J. of Fracture Mech. , 24 (1984) p. 45. 

18. N.I., Muskhelishvili, Some Basic Problems of the Theory of Elasticity. 4th edn., Noordhoff 
International Publishing, Leyden, 1975. 

19. Williams, W.L., ' The stress around a fault or crack in dissimilar media, ' Bull. Seism. Soc. 
Am. , 49, 2 (1959) p. 199. 

20. Bogy, D.B., ' Two edge-bonded elastic wedges of different materials and wedge angles 
under surface traction, " J. of App. Mech. , 38, 2 (1971) p. 377. 


6. References 


138 




21. Wang, S.S. and Choi, l„ ' Boundary layer effects in composite laminates: Part I - Free edge 
Stresses, * J. App. Mech. , 49 (1982) p. 541. 

22. S.G., Lekhnitskii, Theory of Elasticity of an Anisotropic Body, Holden-Day, San Francisco, 
1963. 

23. Wang, S.S. and Choi, I., ' Boundary layer effects in composite laminates: Part II - Free edge 
stresses, " J, Add. Mech. , 49 (1982) p. 549. 

24. Wang, S.S., Elasticity solution for a class of composite laminate problems with stress 
singularity, In Mechanics of Composite Materials , Edited by, Z., Hashin and C.T., Herakovich, 
Pergamon Press, NY, 1982. 

25. Delale, F., ' Stress singularity in bonded anisotropic materials," Int. J. Solids Structures , 
20, 1 (1984) p. 31. 

26. Ting, T.C.T., and Chou S.C., " Edge singularities in anisotropic composites, " Int. J. Solids 
Structures , 17, 1 (1981) p. 1057. 

27. Bogy, D.B., * The plane solution for anisotropic elastic wedge under normal and shear 
loading, " J. of Add. Mech. , 39, 4 (1972) p. 1103. 

28. Carpenter, W.C., " A collocation procedure for determining fracture mechanics parameters 
at a corners, " Int. J. of Fracture , 24 (1984) p. 255. 

29. L.E., Hulbert, The Numerical Solution of Two-Dimensional Problems of the Theory of 
Elasticity. Ohio State University, Eng. Exper. Station, Bull. 198, Columbus Ohio. 

30. Carpenter, W.C., " The eigenvector solution for general corner of finite opening crack with 
further studies on the collocation procedure, * Int. J. of Fracture , 27 (1985) p. 36. 

31. Sinclair, B.G., Okajima, M. and Griffin, J.H., * Path independent integrals for computing 
stress intensity factors at sharp notches in elastic plate, " Int. J. of Num. Meth. in Engg. , 20 
(1984) p. 999. 

32. Wang, S.S. and Yuan, F.G., " A hybrid finite element approach to composite laminate 
elasticity problems with singularities, " J. of App. Mech. , 50, 4 (1983) p. 1. 

33. Jones, R. and Callinan, R.J., " On the use of special crack tip elements in cracked elastic 
sheets, " Int. J. of Fracture , 13 (1977) p. 51. 


6. References 


139 



34. The International Mathematics Subroutine Library, ZANLYT Subroutine, IMSL, INC. 
Houston, Texas. 

35. Muller, D.E., ' Method for solving algebraic equations using an automatic computer, ' 
Math. Tables and Comput. , 10 (1956) p. 208. 

36. The International Mathematics Subroutine Library, ZSCNT Subroutine, IMSL, INC. Houston, 
Texas. 

37. G., Dahlquist and A., Bjorck, Numerical Methods, Prentice-Hall, Englewood Cliffs, New 
Jersey, 1974, pp. 196-201. 

38. J.N., Reddy, An Introduction to the Finite Element Method, McGraw-Hill Book Co., N.Y., 
1984. 

39. Tong, P. and Pian, T.H.H., * The convergence of finite element method in solving linear 
elastic problems, ' Int. J. Solids Structures , 3 (1967) p. 865. 

40. W.D., Whetstone, EAL Engineering Analysis Language Reference Manual, Engineering In- 
formation Systems, Inc., San Jose, California, 1979. 

41. V.V., Novozhilov, Foundations of the Nonlinear Theory of Elasticity. Graylock Press, 
Rochester, N.Y., 1953. 

42. Pian, T.H.H., ' Derivation of element stiffness matrices by assumed stress distributions, " 
AIAA J. , 2, 7 (1964) p. 1333. 

43. Tasi, S.W. Composite Design 1986 . Thinks Composites, 1986, Dayton, Ohio. 

44. Curry, J.M., Johnson, E.R., and Starnes, Jr., J.H., ' Effect of ply drop-off on the strength of 
graphite-epoxy laminate,' Virginia Tech Center for Composite Materials and Structures Report 
CCMS-86-07 and College of Engineering Report VPI-E-86-27, Blacksburg, Dec. 1986. 

45. Whitney, J.M. and Nuismer, R.J., ' Stress fracture criteria for laminated composite con- 
taining stress concentrations,' J. of Composite Materials , 8, (1974), p. 253. 

46. Kin, R.Y. and Soni, R., 'Experimental and analytical studies on the onset of delamination 
in laminated composites,' J. of Composite Materials , 18, (1984), p. 70. 


6. References 


140 



Appendix A. Material Constitutive Relations 


A.1 Transformation Relations 


The geometric relations between the lamina principal 1-2-3 coordinates system and the 
laminate x-y-z system are depicted in Figure A.1. The principal lamina 1-2-3 material coordi- 
nates correspond to the transverse, normal, and longitudinal to the fiber directions respec- 
tively. The angle <p measures the angle between the 3 and z axis, a positive rotation 
corresponding to rotation of the fiber from the z axis to the x axis. The well known lamina 
constitutive relations in the 1-2-3 system are written as follows: 



[A.1] 
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where, 


S„ = , i = 1,2,3 
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and, S () = Sji . In inverted form the lamina constitutive relation in the 1-2-3 system takes the 
following form; 
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C 11 

Cl 2 

C 13 

0 

0 

0 

Cl 2 

c 22 

C 23 
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0 

0 

C 13 

C 23 

C 33 

0 

0 

0 

0 

0 

0 

C 44 
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0 

0 

0 

0 

0 

C 55 

0 

0 

0 

0 

0 

0 

C 66 



[A.2] 


Cn = 


C 22 - 


C 33 “ 


S 22 S 33 

<0 
CM CM 

CO 

r = 

S 12 S 13 

- 

s 23 s„ 

s 


^23 


s 


S11 S33 - 

q2 

b 13 

p = 

Sl2 §23 

— 

S 13 S 22 

s 


^13 


s 


S 11 s 22 - 

q2 

b 12 

p ^ = 

S 13 S 23 

— 

S12 S33 

s 

1 



s 



Q — _ 1 _ 

'-li 


i = 4,5,6 , 
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The transformation of the stresses and the strains from the 1-2-3 coordinate system to the x- 
y-z system leads to; 


a, = [T,] a x , 


e, = [T 2 ]T x . 


The transformation matrices for the stresses, [T,] and [T 2 ], are given by; 


CTi3 = 


m 

0 

n 2 

0 


0 

1 

0 

0 


mn 0 
0 0 


n 

0 

2 

m 

0 

-mn 

0 


0 

0 

0 

m 

0 

-n 


-2mn 

0 

2mn 

0 

(m 2 - n 2 ) 


0 

0 

0 

n 

0 

m 


[A.3.a] 

[A.3.b] 


[A.4.a] 


and, 


[T 2 ] = 


m 

0 

n 2 

0 


0 

1 

0 

0 


2mn 0 
0 0 


m 

0 

-2mn 

0 


0 

0 

0 

m 

0 

-n 


-mn 

0 

mn 

0 

(m 2 - n 2 ) 


0 

0 

0 

n 

0 

m 


[A.4.b] 


Combining eqs. A.1 through eq. A.4, the lamina constitutive relations in the x-y-z coordinates 
system my be written as; 
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[A.5.a] 


E x — [S] O x , 

a x = [C] e x . [A.5.b] 

The coefficients in the compliance matrix, [Si are 

S„ = m*S 11 + m*n 2 (2S 13 + S 55 ) + n 4 Sj, , 

S ^2 = S 21 = m 2 S ^2 4 n 2 S 23 » 

S 13 - S 31 = (m 4 4 n 4 ) S 13 4 m 2 n 2 (S 11 + S 33 - S 55 ) , 

S 18 = S 51 = m 3 n (2 S 13 - 2 S f1 4 S„) 4 m n 3 (2 S 33 - 2 S 13 - S w ) , 

5 22 = S 22 , 

5 23 = = n 2 S 12 4 rn 2 S 23 ■ 

= S 52 = 2 rn n (S 23 S 12 ) 1 

S 33 = n 4 S 31 4 m 2 n 2 (2 S 13 4- S TO ) 4 m 4 S 33 , 

S w - S 83 = m 3 n{2S 33 - 2S 13 - S w ) 4 mn 3 (2S 13 — 2 3,, + S w ) t 
S 44 = m 2 4 n 2 S w , 

S* = S M = m n {S* - S w ) , 

S ra = 2 m 2 n 2 (2 S 33 - 2S„ - 4S 13 - S 55 ) 4 (m 4 + n 4 ) S 8 , 

S w = n 2 S M 4 m 2 S M , 

and all other coefficients are zero. The coefficients in the lamina stiffness matrix, [C] are 
given as; 

C„ = m 4 C u 4 m 2 n 2 (C 13 + 2C*) 4 n 4 ^ , 

C 12 = C 21 = m 2 C 12 4 n 2 Qa , 

C 13 = C 31 - (m 4 4 n 4 ) C 13 + m 2 n 2 (C 11 + C 33 - 4C 58 ) , 

C 13 = C 5 -J = m 3 n (C 3 3 C 33 ) 4 m n 3 (C 33 C 13 ) , 

C 22 = C a , 

C 23 = C 32 = n 2 C 12 4 m 2 C 23 t 

C 2 5 52 C 82 “ m n (C 23 — C 12 ) , 

C 33 = n 4 C t1 4 2 m 2 n 2 (C 13 4 2g 4 m 4 C„ , 
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2 m n (m 2 - n 2 ) C, 


II 

8 

IO 

C 53 — m n 3 (C 13 

Cn) + 

m 3 n (C33 C 13 ) 

c« = 

m 2 C44 4 * n 2 C B8 , 



c« = 

C w = mn(C u - 

c«) . 


C55 = 

2 m 2 n 2 (C 33 + C t1 - 

1 

« 

d 

CNJ 

1 

2 C M ) + (m 4 + 

Ol 

S 

II 

n 2 C44 *4 m 2 C 68 , 




and all other coefficients are zero. 


i 


4 


A.2 Integrated Material Properties 


The integrated laminate properties are obtained by smearing the individual lamina prop- 
erties throughout the thickness of the laminate. This is achieved by defining an average stress 
through the laminate thickness, h, i.e., 



[A.6] 


Substitution of eq. A. 5 .a into eq. A.6 leads to the laminate constitutive relations: 



[A. 7 .a] 


where the overbar is dropped from the stresses for convenience. The laminate stiffness 
components A„ are given by, 
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4 


4 


4 


4 


4 


4 


4 


4 


[A.7.b] 


A ij “ -j“ £ C, i ^ yk ~ ^k-1 ) • 


n being the number of laminae and y k and y k _, are defined in the same manner as in classical 
laminate theory (CLT). In inverted form, 

a^ a 12 3^ 0 0 0 

3i2 322 ®23 0 0 0 

a 13 a 23 a 33 0 0 0 

0 

0 0 0 0 a 33 o 

0 0 0 0 0 a R , 


V 
y xz 

Yxy 


0 0 a u 0 0 


l yz 


l XZ 


[A.8] 


l xy 


It should be noted that although each individual lamina is considered to be anisotropic, the 
laminate constitutive law is that of orthotropic material. The coefficients of the laminate com- 
pliance matrix , a IJt are given by; 


a ii - 


a 33 “ 


II 

> 

IO 

to 

A 33 

CO 
CM CM 

< 

1 

CO 

CM 

03 

CO 

sr* 

< 

CM 

< 

II 

- 

A 23 A 11 


A 

t 


A 


_ A 11 

A 33 

“ A 13 

a 13 

_ A 12 A 23 

— 

A 13 A 22 


A 

t 


A 


_ A 11 

A 22 

.2 

- a 12 

a 12 

_ A 13 A 23 

— 

A 12 A 33 


A 

I 


A 



a li 

A ii 

, i 

= 4,5,6 , 




were, 


A — A^ A 2 2 A 33 


An A 23 A 22 A 13 “ A 33 A 12 + 2 A 12 A 13 A 23 
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Appendix B. The Eigenvalue Problem 


The eigenvalue problem for X and 8 is developed in this appendix for a general bimaterial 
composite wedge depicted by Figure B.l. The traction-free wedge boundaries are located by 
the angles a, and a* in material 1 and 2, respectively. 


B.l The 5 Eigenvalue Problem 


The set of equations which form the eigenvalue problem for 5 is written symbolically as, 

= 0 , [B.l] 

where the coefficients q, } are given by, 

q„ = (cos a, + i u (1) sina 1 ) 4+1 , q,* = (cos a, — iu (1) sina 1 ) 4+1 , 

dl3 = dl4 = d21 = d22 = 0 • 
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q» = ( cos a 2 - io® sina 2 )* +1 , q 24 = (cosa 2 + I u® sina 2 ) s+ ’ , 

^31 = ^32 = 1 • ^33 = P34 = ~ 1 . 

<Ul = - q 4 2 = ^ (1) . $43 ■ - P44 - V 21 • 


with, 




fig. 

m 


k = 1,2 , 


and. 


v R| - 1^1 


t(k) R (k) 


P& . k = 1,2 


The eigenvector , D , is given by, 


D = 


)(D 

;<2> 


[B.2] 


where the two eigenvectors, D® and 6® are related to material 1 and 2, respectively. Due to 
the boundary conditions, the two vectors are related by an equation of the form 





[B.3.a] 


leading to, 


D< 2 > = 


2q 


( 2 ) 


d® + n<”) n® - V”)‘ 
01® - n (,| i <1® + V 1 ') 




The reduced eigenvalue problem takes the form, 
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fl 


4 


4 


I 


4 


4 


4 


4 


4 


4 



where, 


jn = ( cos a, + i u (1) sin a 1 ) 5+1 , j 12 = ( cos a, - i u (1) sin a 1 )‘ +1 , 

j 2 i ~ TrV t (11® + q w )( cos a 2 - i u® sina 2 )* +1 + (q® - r| (1) )( cos a 2 + i u® sin a 2 )* +1 ] , 

2 t|' z ' 

in - [ (q® ~ q®)( cos a 2 - i u® sina 2 ) s+1 + (q® + q®)( cos a 2 + i u® sin a 2 )‘ +1 ] , 


B.2 The X Eigenvalue Problem 


The eigenvalue problem associated with X is written symbolically as, 

[ A {a 1t a 2 ; A.) ] C = 0 , [B.5] 

where, 

a 1k = ( cos a, + pj, 1 ’ sin a 1 ) u+2 > , a 1(k+4) = 0 , 

a*k = (pi 11 cosa, - sin a 1 )( cos a, + p(?» sin a,) tt+1 » , a 2(k+4) = 0 , 

a 3 k = 0 , a 3(k+4) = (cosctj - p|?> sin a 2 ) (i+a , 

a* = 0 , a 4(k+) = (pjficosctj + sinaj)( cos a 2 - pj?> sin a 2 ) (x+1 > , 

a 5 k = 1 , a 5(k+4 , = —1 , 
a «k = Up’ , a s(k+4) = - p|*> , 

a 7 k = Pi 1 ’ . a 7(k+4) = - pj« , 

a«k = qf* . a, (k+4l = - qjf> , 

where k = 1,2, 3, 4, and p k and q k are defined by 

Pk = Pi 1 Pk + Pi 2 > 
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In addition, n k has one of four unique value given by the characteristic equation 

P 22 + (2 Pm + Pee) + 0ii Mk — 0 • 

Such roots do exist and they are always complex or imaginary. 

The eigenvector C is define by, 



[B.6] 


where, the C (1) and C <2> are two vectors associated with material 1 and 2, respectively. Due to 
the boundary conditions, the two vectors are related by an equation of the form 

[E (1) ]C (1) = [E (2) ]C (2) , [B.7] 


where, 


e 1k - 1 


JO = JO JO = JO 

e 2k Pk - e 3k Pk 


JO = JO 

e 4k Mk 


1 — 1.2 . k — 1,2, 3, 4 


The relation between the two vectors is then given by, 

c (2 > = [t]c (T) , 


[B.8] 


where, 

[T] = [E (2) r 1 [E (1) ] , 

is the vector transformation matrix. This leads to the reduced eigenvalue problem associated 
with C (1 > only. The problem may be written as 
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• 



• 

[H(a 1l o 2 ;X)]C (1) = 0 , 

[B.9] 


where, 


# 

h,j = a,j , i = 1,2 , j = 1,2, 3,4 , 


• 

4 

fyj = a i(k + 4)*kj . i = 3,4 , j = 1,2, 3, 4 , 



where t*, are the elements of the transformation matrix [T], 


• 

B.3 Eigenvector Representation 


• 

The eigenvectors D n and C„ associated with the n" 1 eigenvalue, given by eqs. 

B.2 and B.6, 


respectively, are of arbitrary magnitude. To eliminate some of the arbitrariness the vectors 

• 

are normalized according to the following procedure. Letting D{J = D n ( — j — , 

V 2 
to, 

—j — ) leads 
V 2 

• 

5<’> . D ( d " 1 ' | 

" n j 43 J ■ 

[B.10] 


where, 


• 

d (1) - 1 + i 1 

V? VF ' 


• 

d «> _ 

“ “ "d7 ' 
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• 





and D„ is an arbitrary constant, D„ = d„ + i d'„ . If 5 is real, then the normalization of the vector 
DJ, 11 is achieved by letting d},V = 1.0 and D„ = d„ . The eigenvector associated with material 
2, D® , is obtained from the vector transformation eq. B.3.b. 

In a similar manner letting, 


leads to, 


where, 


p(t) _ p / 1 

Sil c n I 1 = 


V? ' Vi - ' 



-PI 

C n1 





I = 2,3,4 , 


Cb.ii] 


and C„ is an arbitrary complex constant, C„ = c n + i c'„ . As before, if X is real, then 
c{,V = 1.0 and C„ = c n is a real constant. The eigenvector C® associated with material 2 is 
obtained via the vector transformation eq. B.8. 


The unknown coefficients d$ and cJ,V . • = 2,3,4 are obtained by eqs. B.4 and B.9, re- 
spectively, by eliminating the first row and first column from matrix JJ] and [H], respectively. 
Thus, for both 5 and X there is one complex undetermined coefficient associated with each 
eigenvalue. 
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B.4 Eigenfunction Expansion Representation 


The eigenfunction expansion representation of F (,) (x,y) and X F (#) (x,y) will differ depending 
on whether the 8 and X are real or complex. As noted before, if 5 and X are complex they will 
occur in complex conjugate pairs. In the following section the proper functional form for a real 
and complex eigenvalue is discussed. Only detailed development for F 1 * 1 (x.y) is shown since 
the form for 'P , * ) (x,y) is identical. 

Real Eigenvalue 

For a real eigenvalue , 


X = y (real) 

The stresses associated with F (,) (x.y) can be written as, 

2 _ 

CT i = L (A k G|j{ + A k + 2 G lk ) , i = x,y,xy , [B.12J 

k = l 


where, 


A k = a k + i a' k 
G ik = 9| k + ■ 9'ik 

and G lk is the conjugate of G lk . In addition g lk and g', k are unique functions of r and 9 ., as 
shown next. Considering o y for example, 

ct v - I A k Z k . [B.13] 

k = 

By making the following substitutions; 
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X 


r cos 0 , y 


r sin 9 , 


[B.14.a] 


Z = x + |iy = rpe 


[E,14.b] 


where, 


P = a ± i|3 , 


<j y becomes 


p 2 = { cos 0 + a sin 0) 2 + p 2 sin 2 0 


[B.14.c] 


tan y = 


P sin 9 

cos 0 + a sin 


[B.14.d] 


Oy — A 1 Gy 1 + A 2 Gy2 + A3 Gy.) + A4 Gy 2 


[B.15] 


In the above 


G y i = (r p,) 7 [ cos (y Vi) + i sin (y Vi)] . 


G y2 = (r P 2 ) y [ cos (y y 2 ) + i sin (yy 2 )] . 


As noted earlier, p always occurs in a complex conjugates pairs, that is, 


-> AH + ift , pj - <z 2 + ip 2 


Pa = oi - iPi , p 4 " “ i P2 • 


Similar results can be obtained for the other stress components. Expansion of eq. B.12 leads 


a, = X C(a k + a k+2 )g, k + (a' k+2 - a' k )g' ik ] 

k = 1 
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+ iC(a' k + a'k+2)9jk + ( a k — a k+ 2 )fl , ik] • 

However, since <r, must be a real quantity and g ilt and g' lk are unique functions of r and 0 , 
setting the imaginary part of a, to zero leads to, 

a'k+2 = — a 'k • a k + 2 ~ a k • 

that is A k+J = A k . Next the coefficients are redefined, 

c k = a k + a k+2 = 2a k 

c k+2 = a ’k+2 “ a V = ~ 2 a k • 

It follows that the stresses can be written as, 

a, = £ Cc k Re(G, k ) + c k+2 lm(G ik )] , i = x.y.xy , [B.19] 

k = 1 

with c k , (k = 1,2, 3, 4) being real constants. 

Complex Eigenvalue 

As was previously stated, if X is complex it will occur in complex conjugate pairs, i.e., 

X = y + i<p, I = y — i <p . 

The stresses associated with F w (x,y) , for a pair of complex conjugate eigenvalue, will take 
the following form; 

4 4 

<T| = S A k G lk + L B k G (k , i = x.y.xy , [B.20] 

k=1 k=1 

with A k and B k being arbitrary complex constants, A k = a k + i a' k and B k = b k + i b' k . The 
function G ik associated with X is the conjugate of the function G, k = g ik + i g' (k associated 


Appendix B. The Eigenvalue Problem 


157 


with X . In addition, g ilt and g' ik are unique functions of r and 0, as shown next. Considering a y 
again for X and X in the most general form becomes, 

5y = I A k z[ + £ B k zl , [B.21] 

k=1 k=1 

where for complex eigenvalues, using relations B.14.a and B.14.b, 

Z k = (rp k ) T e - Vk<,, [ cos q k + i sin q k ] , [B.22.a] 

Z k = (rp k ) Y e Vlt<p [ cos q' k + i sin q' k ] , [B.22.b] 

where, 

T| k * YYk + <P ln(rp k ) . 

T|'k = YVk “ <P ln(r p k ) . 

Considering eqs. B.14 and B.17 and substituting eq. B.22 into eq. B.21, o r may be written as, 
C7 y = A 1 f., + A 2 f 2 + A 3 f 3 + A 4 f 4 + A 5 f 5 + A 6 f 6 + A 7 f 7 + A 8 f 8 , [B.23] 

where, f, , i = 1,2 8 are complex functions. Next, it can be shown that f 7 is the conjugate 

of f, (i.e., f 7 = ?,) where, 

ft = (rpi) T e' l ' 1<p {cos [yvi + <P Mrp,)] + i sin + <p Infrp.,)]} 

Similarly, f 8 = f 2 , f 5 = f 3 , and f, = f 4 . Since the numbering of the functions and coefficients is 
arbitrary, eq. B.23 can be written in the form of eq. B.20, where, 

Gyl = h < ®y2 = ^2 ■ ®y3 = ^3 • ®y4 = h 

Gyl = fy , = f 2 , G y3 = f 3 , G y4 = f 4 . 
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Similar results can be obtained for the other stress components. Expansion of eq. B.20 leads 
to the following result; 


4 

o, = E C(a k + b k )g lk + (b' k - a' k )g' ik ] 

k = 1 

+ i C (a' k + b' k ) g jk + (a k - b k )g' lk ] . [B.24] 

However, since a, must be a real quantity and gi k and g' lk are unique functions of r and 6. 
Setting the imaginary part of a, to zero leads to, 

b k = — a k • b k = a k • 

i.e., B k = A k . Renaming the constants as, 


c k a k b k 2 b k , 

c 'k = b 'k ” a 'k = 2 b' k , 

leads to the following eigenfunction expansion of stress, 

44 r 
o, = E c k Re (G ik ) + E c' k Im (G jk ) , , i = x.y.xy , LB.25J 

k=l k=1 

where c k and c' k are two unknown real constants. It should be noted that the above repre- 
sentation of stresses leads to the classical stress formulation arising from a complex stress 
function. For example, 


o x = 2 Re [ E C k ] , etc... 

k = i dy 2 

with C k being an arbitrary complex constant. 
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Based on the above discussion and the one given in the previous section, the general form 


of the stress functions F w (x,y) and 4 /(,) (x,y) can be written. 


Real Eigenvalue 


For real n^ eigenvalue, the two stress functions F (,) (x,y) and 'P w (x,y) are given by, 


F w (x,y) = E 


"-1 (X. + 1)(X B + 2) 


{ReCE c nk Z^ + *> ] 


+ lm[ c n(k+2) zl : Xn+2) ] > , 


[B.26] 


^•>(x,y) = £ — { Re [ d n1 Z( 5 n +1) ] 

n- 1 (8 n + 1) 

+ lm[d n2 z{f B+1) ] > , [B.27] 

with c„ and d B being arbitrary real constants and c Bk , (k = 1,2, 3, 4) and d nk , (k = 1,2) known 
quantities of the normalized n ,h eigenvector (see section B.2). 


Complex Eigenvalues 


If the n th eigenvalue is complex, the two stress functions are given by, 


00 4 7jl n + 2) 

F<*> (x,y) = E { c n Re [ E c Bk , 

n-1 k-1 (X n + ^(K n + 2) 


4 Z [X, + !) 

+ c n Im [ E c nk - 


k = i nK (X n + 1)(X n + 2) 


]> . 


[B.28] 


¥<•> (x,y) 


£ { d B Re [ E d 

n-1 k-1 


nk 


Zjfn +2 > 

(5 n + D 


] 
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[B.29] 


2 z<S„+D 

+ d' n lm [ X d nk - — 3 > 

k = i (b n + 1) 

with c n ,c'„,d n and d' n being arbitrary real constants and c nk , (k = 1,2, 3, 4) and 
d nk , (k — 1,2) being known quantities of the normalized n ,h eigenvectors. With the functional 
form for F w (x,y) and T (,) (x,y) known, the eigenfunction expansion for both stresses and dis- 
placements can be determined. 
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Appendix C. Finite-Element Formulation 


C.1 PE2D Finite-Element Formulation 


In this appendix the PE2D finite-element program is described in greater detail. This pro- 
gram is based on the elasticity equations which govern an elastic body for which the stresses 
do not vary along it's generator (the z axis). Such a condition is referred to as the a general- 
ized plane deformation. 

The finite-element program is based on the displacement formulation approach and it is 
derived from the total potential energy function, n p , 

n p = \ j v II 5 X dv - J A u T t dA , [c.i] 

The first term corresponds to the body's total strain energy, the second to the work done by 
the external force, t. The work done by the body forces is neglected. The generalized plane 
deformation finite-element model is based on the displacement field for homogeneous 
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anisotropic body for which the stress do not vary with z for this situation. The displacement 
functions u, v, and w have the following form: 


u = - Bl ^ 3 - z 2 - B 4 y z + U(x,y) + <d 2 z - oo 3 y + u 0 , [C.2.a] 

V _ _ B 2 ^33 _ z 2 + b 4 X Z + V(X,y) + O 3 X - <»! z + v 0 , [C.2.b] 


w = ( B, x + Bj y + B 3 ) a 33 z + W(x,y) 

+ CD, y - co 2 x + w 0 , [C.2.c] 

where B, , I = 1,2,3, 4 , are arbitrary constants of integration, co, , i = 1,2,3 are rigid body ro- 
tations, and u„ , v 0 and w e are rigid body translations. In addition, U(x,y), V(x,y) and W(x,y) are 
unknown displacement functions which depend on x and y only. This displacement field can 
be written in vectorial form as, 


u = u 0 + U 


[C.3] 


where, 


Uq = {u 0 (x,y,z) ; v 0 (x,y,z) ; w 0 (x,y,z)> ; , 
U T = (U(x,y) ; V(x,y) ; W(x,y)> , 


and, 


u 0 (x,y,z) = - Bi 2 333 z 2 - B 4 yz , 

v 0 (x,y,z) = - B2 - 333 z 2 + B 4 xz , 
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w 0 (x,y,z) = (B^ + B 2 y + B 3 )a 33 z 


4 


The terms associated with the rigid body rotation and translation are omitted. It may be shown 
that the unknown constants B, , i = 1,2, 3, 4 , are related to the body's kinematic end condi- 
tions, i.e., 


B 3 = 


_£o_ 

a 33 


b 2 


Jk. 

a 33 


b 4 



In the stiffened-skin structural context, e 0 and k, are the axial extension and curvature in the 
z direction, and k„ is the twist curvature about the z axis. The coefficient B, relates to the in- 
plane twist about the y-axis and is of no consequence in the structure considered here and is 
therefore set to zero. The strain vector, e x is given by, 


e x - E + e 0 , 


[C.4.a] 


where, 


e T = 


{e x I Ey , £ z ; Yyz • Yxz • Yxy) 


[C.4.b] 


p t = r_aj_ . av . n . jw ■ M. ■ iy. + jy_\ 
i dx'dy , ’dy , dx’dy 3x ' 


[C.4.c] 


and, 


{o > o i (©o y) ’ 



-f-y ; o} 


Next, the laminate constitutive law may be written in a reduced form as, 


[C.4.d] 


4 


4 


4 


i 


i 


4 


4 


i 


4 
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4 





Ey = 


dU 

OX 

av 

A 

aw 

ay 

aw 

ax 


- Ey 


e„ 


^yz 

y°xz 


au f av 
ay ax 


Pi 1 Pi 2 0 

Pi 2 022 0 


0 0 
0 0 
0 0 p w 0 0 

0 0 0 p 55 0 

0 0 0 0 Pe6 


u y 

V 

Yxz 

Yxy 


[C.5] 


where, 


E ° = ‘& (e ° + Kzy) : e ° = -& (e ° + KzV) 


and, 


y°z = 


_ K xz 


Y « = 


v xz 


Plj - a lj - 


3ii 3 


13 d j3 
a 33 


i.j = 1 , 2 , 4, 5, 6 . 


Notice that E and e 0 are modified by condensing Oj out (i.e., the reduced compliance form). In 
inverted form eq. C.5 will be referred to by, 


a x = [Q] e x , 


[C.6] 


where, [Q] = [p] -1 

Considering eqs. C.1 and C.5 it may be concluded that the problem associated with U(x,y) 
and V(x,y) is decoupled from the problem which involves W(x,y). In addition, since U, V, and 
W are functions which depend on x and y alone, only a two-dimensional finite-element model 
is required. Considering eqs. C.3, C.5 and C.6 the total potential energy for a typical element 
of thickness t* in the z direction becomes, 
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dA + 




- -rI»- gTQE<IA - '°J s .S T Q^dA + 4 j a . 


( e o + K z y) z 
a 33 


■f &.?QS,dA - t«; s . u T . ds - 


dS 


[C.7] 


where A* is the element surface area in the x-y plane and SJ is the part of the element 
boundary on which external traction forces are acting. In the displacement formulation 
finite-element program, the unknown displacement functions, U(x,y), V(x,y), and W(x,y) are 
approximated using a Lagrangian interpolation functions, 

U(x,y) = I U| Vi(x.y) , [C.8.a] 

i = i 


V(x,y) 


n 


£ vi Vi(x.y) . 

1=1 


[c.8.b] 


W(x,y) 


n 

£ w,V|(x,y) 


1 = 1 


[C.8.c] 


where, u„ v, and w, are the generalized nodal displacements and y, are the two-dimensional 
Lagrangian interpolation function. The interpolation functions \y, assume different forms, de- 
pending on the type of element used. In the present investigation two types of element were 
used; the 8-node rectangular and the 6-node triangular isoparametric elements. The interpo- 
lation functions for these elements have the following form; 


I 


4 


i 


i 


< 


< 


< 


8-Node Isoparametric Element Interpolation Functions 

ViM = - -J-d -5)0 -t|)(§ + ti + i) , 

= j(i + $)d -tOfc-n - D . 
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i 


i 


i 



• 

Vs&*l) = 

-J-0 + 5)0 + n)(& + if - i) . 

[C.9.c] 


VK 4 (tTl) = 

-J-d “ §)(1 + r\)(r\ - %- 1) , 

[C.9.d] 

• 

vste-ii) = 

y(1-* 2 )(1-T|) , 

[C.9.e] 

• 

vs(i'n) = 

y(1-H 2 )(1 + §) ■ 

[C.9.f] 



y<1 “ & 2 )d + il) . 

[C.9.g] 

• 

vsU.'n) ■ 

yd -n 2 )d -§) . 

[C.9.h] 



[C.IO.a] 

[C.IO.b] 

[C.IO.c] 

[C.IO.d] 

[C.IO.e] 

Cc.io.f3 


where the elements nodal numbering and natural coordinates system are shown in Figure C.l 
Substitution of eqs. C.8 into eq. C.4.b leads to 
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E = Bq , 


Cc.li] 


where, 


q = [u, V| W|] , i = 1,2 n 


are the nodal displacements. The matrix B is given by. 


B = [B 1 B 2 ... B n ] , 


[C.12] 


where, 


B, = 


0 

ax . 

0 *1 

ay 

0 0 


, i = 1,2 n 


avi _£vi_ 

ay dx 


Finally, substituting eqs. C.2, C.8, and C.9 into eq. C.7 and minimizing the total potential energy 
for a typical element, 


- [t e J A . B T QB dA] q - t e J A . B T Qe 0 


*© r itJ 


k r 


t dS = 0 


[C.13] 


Next if only the solution for U and V are considered, the above integral equations may be cast 
into a set of simultaneous equations in a matrix form 


k 11 k 12 


k 21 k 22 
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or in the more common form, 


[k]q = F . [C.14.b] 

The components of the stiffness matrix, [k] and force vector, F are given as follows: 




3yj 


dy 3y 


[C.15.a] 


u 12 _ ^21 _ + e c T 5 V| ^ , 

k'j - klj - t Ja« [ Q 12 -^ + Qes-^-^Jdxdy . 


[C.15.b] 


* - ■ 


[C.15.c] 


F) 1 - t a J A « 0 a 33 ~ ^ a i3 + Qi 2 a 23 ] dx dy 


[C.15.d] 


+ t 6 J s . Vi t x ds , 


F ? - tC J A * a 33 Z ^ ^ Q 12 a i3 + Q 22 a 23 ^ dx dy 


[C.15.e] 


+ t6 Js» V|ty ds , i.j — 1.2 n 


Using the interpolation functions, eqs. C.9 and C.10, the elements stiffness matrices and force 
vectors in eqs. C.15.a through C.15.e can be evaluated numerically. Once [k] and F are known 
the solution for q can be obtained. 
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C.2 EAL Finite-Element Program 


The EAL finite-element program is based on hybrid element formulation. In the displacement 
formulation finite-element approach, the element stiffness matrices are derived using the 
principle of minimum potential energy. Whereas, in the hybrid finite-element formulation the 
element stiffness matrix are derived using the complementary energy principle given by, 

n c = yJ v a T CadV - J A u T t dA , [C.16] 

where the work done by the body forces were neglected. In eq. C.16, C is the material stiffness 
matrix, u is the displacement vector, and t is the surface traction vector. In the hybrid element 
formulation the stresses are written in terms of generalized stress coefficients, 

o = P b , [C.17] 


where, 

b T = [b 1 b 2 ... b n ] , 

and P is a matrix which depends on position and is chosen in such a way so that a satisfies 
the stress equilibrium equations identically. Using eq. C.17 the body internal strain energy 
becomes, 


= y b T {J v pT CP dv } b . [C.18] 

Next the displacements are written in terms of the element generalized nodal displacements, 

u = Lq . [C.19] 
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where L is a function of position and insures inter-element compatibility. In addition, the sur- 
face tractions may be related to the body stresses by, 

t - Rb , [C.20] 

where R is a function of the boundary position. Upon substitution of eqs. C.18 through C.20 into 
eq. C.16 the total complementary energy becomes, 

n c = jb T Hb - b T Tq , [C.21] 


where, 


H = J v P T C P dV , 

T = J A R T LdA . 

The next step involves the minimization of n e with respect to the generalized stress coeffi 
dents, b„ 


dn c 

<3b| 


0 


leading to a relation between b and q, 

b = H -1 Tq 

From which the element stiffness matrix is calculated, 

[kj = T t H - 1 T . 

The solution to the finite-element problem is obtained in the standard form, i.e., 

[k] q = F , 
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except that [k] is replaced by [kj. If [k] is the stiffness matrix of perfectly conforming ele- 
ments, than [kj approaches [k] from the soft side. Whereas, [k], as derived from the dis- 
placement formulation finite-element, approaches [k] from the stiff side. 


EAL - E43 and E33 Elements 


The E43 and E33 are the 4-node rectangular and 3-node triangular plate elements. Each 
node of these elements has 5 degrees of freedom (dof) consisting of two rotation and three 
displacement components. The element formulation is based on the assumed stress resultant 
polynomials which satisfy the following plate equilibrium equations; 


fQx 

8x 


+ 



q 


o , 


dx 


3y 


Qx 


dx dy 


?*** + 

sx 2 sxay ay 2 5 


0 , 


where q is applied uniform transverse load. All symbol and coordinate norrffenclature are 
those common to the theory of plates. The stress resultant polynomials which satisfy the 
above equations have the following form; 
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E33 Element 


N x = b 3 


Ny b 3 , N X y 


bi 




E44 Element 



N X y 


* 

where, 


b 17 



The plate constitutive relations for composite materials are given by, 
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EAL - S61 and S81 Solid Elements 


The S61 and S81 are 6-nodes, pentahedron, and 8-nodes, hexahedron, solid elements. Each 
node of these elements has 3 degrees of freedom, u, v, and w. For these elements the as- 
sumed stress polynomials have the following form; 

a x = b, + b 7 y + b,z + b 1g yz , 

a y = b 2 + b,x + b 10 z + b 17 xz , 

o t = b, + b„x + b 12 y + b 18 xy , 

V = b 5 + b 14 x , 

= b s + b 15 y , 

= b 4 + b 13 z , 

These set of assumed stress polynomials satisfy the stress equilibrium equations three- 
dimensional elastic body. The element constitutive law assumes the most general anisotropic 
form, 


e x [S]a x , 


where [S] is a fully populated matrix. 
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Appendix D. Material Properties and Eigenvalue 
Data 


The material properties used in the current investigation are representative of materials 
used by the aircraft industry. Though not representing any particular material, they are used 
here to demonstrate the analytical model. The lamina principal material properties are, 


E 1 = E 2 = 2.1 Msi , E 3 = 20.0 Msi , 

G 31 = G 32 “ ^12 ~ 0.85 Msi , 

V 31 ■ v 32 = V 12 = 0.21 , 

with the 1, 2, and 3 being the transverse, thickness, and fiber-direction respectively. The inte- 
grated laminate material properties are obtained by smearing the individual lamina properties 
through the thickness of the laminate. Table D.1 lists the engineering constants for the differ- 
ent laminates used in the investigation. In the study, two sets of skin-stiffener material system 
combinations were used. In both of these combinations the skin laminate was always main- 
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Table D. 1. Laminate Material Properties. 





























































tained at ( ± 45/90*),. The stiffener flange laminate layup varied. The following four stiffener 
flange laminates were used: 1. ( ± 45/0/90), ; 2. ( ± 45/0/90/0/ ± 45), ; 3. ( ± 45/0*), ; and 4. 
( ± 45/0/ ± 45),. In the stiffener parametric study context, laminate 1 is the baseline stiffener 
flange configuration, laminate 2 is the 'thick flange", laminate 3 is the 'soft flange', and lami- 
nate 4 is the 'thick & soft flange" configuration. From Table D.1 it is clear that the combination 
of the skin laminate with either the first or the second flange laminate (i.e., laminate 1 or 2) 
will produce the same set of eigenvalues for a given flange termination angle. This material 
system combination is referred to as combination no. 1. The same is true for skin-stiffener 
material combination which use flange laminates 3 or 4. This material system combination is 
referred to as combination no. 2. The eigenvalues for these two material system combinations 
and for three flange termination angles are given in Table D.2 and D.3, respectively. 

Finally, the adhesive layer material properties used to produce the results in section 2.6 
were taken from Delale and Egdogan 1 and are listed below; 

E = 0.445 Msi , 

G = 0.165 Msi . 


1 Delale, F., Erdogan, F., and Aydinoglu, M.N., ' Stresses in adhesively bonded joints: A closed-form 
solution,' J. of Composite Materials , 15, (1981) p. 249. 
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Table D.3. The First 20 Eigenvalues for Material System Combination no. 
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